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We discuss a certain class of two-dimensional quantum systems which exhibit conventional 
order and topological order, as well as two-dimensional quantum critical points separating 
these phases. All of the ground-state equal-time correlators of these theories are equal to 
correlation functions of a local two-dimensional classical model. The critical points therefore 
exhibit a time-independent form of conformal invariance. These theories characterize the 
universality classes of two-dimensional quantum dimer models and of quantum generaliza- 
tions of the eight-vertex model, as well as Z2 and non-abelian gauge theories. The conformal 
quantum critical points are relatives of the Lifshitz points of three-dimensional anisotropic 
classical systems such as smectic liquid crystals. In particular, the ground-state wave func- 
tional of these quantum Lifshitz points is just the statistical (Gibbs) weight of the ordinary 
2D free boson, the 2D Gaussian model. The full phase diagram for the quantum eight-vertex 
model exhibits quantum critical lines with continuously- varying critical exponents separating 
phases with long-range order from a Z2 deconfined topologically-ordered liquid phase. We 
show how similar ideas also apply to a well-known field theory with non-Abelian symmetry, 
the strong-coupling limit of 2 -|- 1-dimensional Yang-Mills gauge theory with a Chern-Simons 
term. The ground state of this theory is relevant for recent theories of topological quantum 
computation. 

I. INTRODUCTION 

During the past decade and a half there has been an intense search for new kinds of theories 
describing quantum condensed-matter systems. Many experimental results have implied that strongly- 
correlated fermionic systems exhibit qualitatively new types of physical behavior. The now-classic 
example of this is the fractional quantum Hall effect, where one of the striking consequences of strong 
correlations is that the Laughlin quasiparticles have fractional charge and fractional statistics, even 
though though the microscopic degrees of freedom are electrons with integer charge and fermionic 
statistics 

Traditionally one classifies different phases in terms of order parameters which give a global char- 
acterization of the physical state. In turn, the local fluctuations of this order parameter field drive the 
phase transitions between ordered and disordered states of these systems. This viewpoint, pioneered 
by Landau and his school, has been extremely successful in condensed matter physics and in other 
areas of physics, such as particle physics, through the powerful underlying concept of spontaneous 
symmetry breaking. Much of the structure of modern theory of critical phenomena is based on this 
point of view 

However, there are many different experimentally-realizable phases (and even more realizable the- 
oretically!) in the fractional quantum Hall effect, but no local order parameter distinguishes between 
them. These phases are incompressible liquid states which have a fully gapped spectrum and do not 
break any symmetries of the Hamiltonian. The lack of a local order parameter led to many interesting 
discussions of the off-diagonal long-range order in the Hall effect One particularly elegant way of 
characterizing the order in the fractional quantum Hall effect is as topological order [5|. The topologi- 
cal order parameters are non-local; they are expectation values of operators which are lines or loops. 
Because of this, they can (and do) depend on topology: their value depends on the genus of the two- 
dimensional surface on which the electrons live. One interesting characteristic of a topological phase is 
that the correlation functions in the ground state do not depend on the locations of the operators, but 
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only on how the loops braid through each other. In addition, the degeneracies of these non-symmetry 
breaking ground states on topologically non-trivial manifolds are determined by the topology of these 
manifolds 6] . 

Although so far the only unambiguous experimental realizations of topological phases are in the 
fractional quantum Hall effect, there has been considerable effort to find, both theoretically and ex- 
perimentally, condensed matter systems whose phase diagrams may exhibit topological ground states. 
Much of the current work involves studying fractionalized phases in time-reversal invariant systems 
(see e.g. 0@|)- One reason is that the "normal state" of high-temperature superconductors lacks an 
electron-like quasiparticle state in its spectrum. There are reasons to believe that frustrated magnets 
may also exhibit fractionalized behavior as well. 

A particularly well-known and simple model with a topological phase is the quantum dimer model, 
which was invented as a way of modeling the short-range resonating-valence-bond theory of super- 
conductivity 0]. The degrees of freedom of this two-dimensional model are classical dimers living 
on a two-dimensional lattice. With a special choice of Hamiltonian (called the RK point), the exact 
ground-state wave function can be found When the dimers are on the square lattice, the result 
is a critical point. If one deforms this special Hamiltonian, one generically obtains ordered phases. 
However, a topological phase occurs in the quantum dimer model on the triangular lattice [^. When 
the quantum dimer model is in a topological phase, an effect analogous to fractionalization occurs 0. 
This is called spin-charge separation. One can view the dimers as being created by nearest-neighbor 
pairs of lattice electrons in a spin-singlet state. Even though the fundamental degrees of freedom (the 
electrons) have both spin and charge, one finds that the basic excitations have either charge (holons) or 
spin (spinous), but not both. To prove this occurs, one must show that if one breaks apart an electron 
pair (dimer) into two holons or two spinous, they are deconfined. For the triangular-lattice quantum 
dimer model, this was shown in Ref. 8]; the analogous statement in terms of holon-holon correlators 
was proven m Ref. m. At a quantum critical ph^se transition between an ordered/confining phase 
and a disordered /deconfini ng pha se (or between different confining states), confinement is lost: the 
RK point is deconfining jl2lll3ll4l |. 

The notion of spin-charge separation is one of the basic assumptions behind the RVB theories of 
high-temperature superconductivity [9l ll5lil^il^il^il9ll20ll2lj|22| . which effectively can be regarded as 
strongly-coupled lattice gauge theories. In 2 + 1-dimensional systems spin-charge separation can only 
take place if these gauge theories are in a deconfined phase |il3ii23>»24] . In 2 + 1-dimensions this is 
only possible for discrete gauge symmetries. For a continuous gauge group, say U{1) or SU{2), 2 + 1- 
dimensional gauge theories are always in a confining phase, unless the matter fields carry a charge 
higher than the fundamental charge so that the gauge symmetry is broken to a discrete subgroup 
[25,,] . Thus, the only consistent scenarios for spin-charge separation necessarily involve an effective 
discrete gauge symmetry, which in practice reduces to the simplest case Z2. Of particular interest is 
the fact that the low-energy sector of the deconfined phases of discrete gauge theories are the simplest 
topological field theories 26 . 27,2^] . 

Many of these ideas have their origin in the conceptual description of confined phases of gauge 
theories as monopole condensates, and of their deconfined states as "string condensates" In gauge 
theories it has long been known that their phases cannot characterized by a local order parameter, 
since local symmetries cannot be spontaneously broken. The phases of gauge theories are understood 
instead in terms of the behavior of generally non-local operators such as Wilson loops and disorder 
operators [sollsH . a concept borrowed from the theory of the two-dimensional Ising magnet [s^l- 

Interesting as they are, the applicability of these ideas to the problem of high-temperature super- 
conductivity and other strongly correlated systems is still very much an open problem. Topological 
fractionalized ground states are not the only possible explanation of the unusual physics of the cuprates. 
In fact, when constructing local microscopic models of strongly-correlated systems which are suspected 
to have fractionalized phases, many theorists have found that instead these models have a strong ten- 
dency to exhibit spatially-ordered states, a.k.a. "valence bond crystals", which appear to compete 
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with possible deconfined states. It is now clear that the regimes of strongly-correlated systems which 
may favor fractionalized phases also favor, and perhaps more strongly, non-magnetic spatially ordered 
states of different types, including staggered flux states or d-density wave states js^], and elec- 
tronic analogs of liquid crystalline phases 0111. By now there are a number of examples of models 
with short-range interactions whose phase diagrams contain both fractionalized and spatially-ordered 
phases {8i>36i.37.i38il] . It has recently been proposed that deconfined critical points may describe the 
quantum phase transitions between ordered Neel states and valence bond crystals [s^ . 

For several reasons, most of the studies of topological order and quantum critical points have 
focused on examples with two spatial dimensions. The experimental reason is that the Hall effect 
is two-dimensional, and typical strongly-correlated systems, such as the cuprate high-temperature 
superconductors, are often effectively two-dimensional. Theoretically, it is because in two dimensions 
particles can have exotic statistics interpolating between bosonic and fermionic A common 

characteristic of topological phases in two dimensions is the presence of exotic statistics, which occur 
in the fractional quantum Hall effect _4j]. The statistics can even be non-abelian: in some cases, the 
change in the wave function depends on the order in which particles are exchanged Systems 
with non-abelian statistics are particularly interesting because they are useful for error correction in 
quantum computers ^42,.43..44.45.46] . 

In this paper we will discuss models with topological phases and ordered phases, as well as quantum 
phase transitions separating them. There has also been a great deal of interest in quantum critical 
points in and of themselves [43]. At a quantum critical point, the physics is of course scale invariant, 
but it need not be Lorentz invariant. The quantum critical points discussed in this paper have dy- 
namical critical exponent z = 2, instead of the usual z = 1 of a Lorentz-invariant theory. This allows 
for some striking new physics. The action of these z = 2 quantum critical points is invariant under 
time-independent conformal transformations of the two-dimensional space. A remarkable consequence 
is that the ground-state wave Junctionals of the field theories discussed here are conformally invariant 
in space. This means that the ground state wave functional is invariant under any angle-preserving 
coordinate transformations of space. For two-dimensional space, there is an infinite set of such trans- 
formations, as is familiar from studies of two-dimensional conformal field theory |4^. This sort of 
behavior is not common at all: the action of a field theory at a critical point is often scale invariant 
(and also conformally invariant), but the ground-state wave functional itself in general is not. We dub 
critical points with this behavior conformal quantum critical points. 

One of the consequences of the conformal invariance of the ground state wave function is that all the 
equal-time correlators of the quantum theory are equal to suitable correlation functions of observables 
of a two-dimensional Euclidean conformal field theory. We will exploit this connection in this paper 
quite extensively. However, just as important, conformal invariance of the wave function implies that 
the ground state of this 2 -|- 1-dimensional theory at a conformal quantum critical point must have 
zero resistance to shear stress in the two-dimensional plane. 

This can be seen as follows. Consider an infinitesimal local distortion of the geometry of the 
two-dimensional plane represented by an infinitesimal change 5gij{x)ol the two-dimensional metric, 
as is conventional in the theory of elasticity Let |^) be the ground state wave function for the 
undistorted plane and \^{g)) be the ground state wave function in the distorted plane with two- 
dimensional metric gij{x) = 5ij-\-5gij{x). Under this distortion the Hamiltonian of the system changes 
by an amount 



6H{g) = I d'x j—^ 6gij{x) + ... (1.1) 



^ Recall that the change in the metric, given by the strain tensor, is quadratic in the local deformation of the system. 
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To first order in perturbation theory in 5H, the change of the ground state wave energy is 

where Eq is the exact ground state energy of the distorted system. On the other hand, the chan ge o f 
the norm of the ground state wave function ||^'|| is, to all orders in perturbation theory, given by|49| 

m' = ^-^ (1.3) 

where eo is the ground state energy of the undistorted system. Thus, the change of the norm ||^'|| is 
determined by the (2 x 2) stress tensor Tij of the 2 + 1-dimensional theory 



5H 

On the other hand, we can regard as the partition function Z of a two-dimensional Euclidean 

pcfl 
ij 



T^3i^) = {T:rTzx) ■ (1-4) 



conformal field theory. This theory has an Euclidean stress-energy tensor, T^J^, defined by[4i 



^* = -T^ =f^> (1-5) 



\Sgij{x)/ Sgij{x) 

which essentially coincides with the stress tensor of the 2 + 1-dimensional quantum field theory de- 
fined above. Scale invariance, rotational invariance and conservation require that Tij be a conserved 
(divergence free) symmetric traceless tensor. Consequently, the effective Hamiltonian (as well as the 
action) at this quantum critical point can depend on the spatial gradients of the field only through 
the "spatial curvature", e.g. (V'^(p)'^ in a scalar field theory. In other words, at a conformal quantum 
critical point for a scalar theory, the stiffness vanishes: the usual (Vt/?)^ term is not possible. This 
means that the dynamical critical exponent of this quantum critical theory must be z = 2. We call 
such theories quantum Lifshitz theories; we will discuss such critical points in detail. 

In this paper we discuss lattice models which exhibit both ordered/confined phases and disor- 
dered/deconfined phases. We will also discuss the field-theory description of these phases and of the 
phase transitions. To simplify matters, and to be able to obtain exact results, we will introduce mod- 
els whose ground-state wave function will be known exactly and whose properties we will be able 
to determine quite explicitly. In this sense, these models are a generalization of the quantum dimer 
model at the RK point. The basis of the Hilbert space of these models is the configuration space of a 
two-dimensional classical statistical-mechanical system or Euclidean field theory. Each of these basis 
states is defined to be orthogonal with respect to the others. An arbitrary state in this Hilbert space 
can therefore be described as some linear combination of these basis elements. Describing the Hilbert 
space in such a fashion is not particularly novel. The unusual feature of the models we will discuss is 
that the ground-state wave function can be expressed as in terms of the action or Boltzmann weights 
of a local two-dimensional classical theory. The normalization of the wave function will then be the 
partition function or functional integral of the classical two-dimensional model. This special property 
is why the wave functionals at the critical points are have a time-independent conformal invariance at 
their critical points. The field theory of these conformal quantum critical points can be extended to 
describe nearby ordered and disordered phases, including their confinement properties. We will study 
this quite explicitly in a quantum generalization of the eight-vertex model. However, much of the 
physics we discuss should apply to topological phases and {z = 2) quantum critical points in general. 

We will also study theories with a continuous non-abelian symmetry. We show that, interestingly 
enough, it is very difficult to construct a non-trivial conformal quantum critical point with such a 
symmetry. We do find a Hamiltonian whose ground state is the doubled Chern-Simons theory of Ref. 
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|50||5l|. This is a time-reversal invariant theory of interest in topological quantum computation and 
in (ordinary) supercondutivity; it is in a gapped topological phase. 

In section^we discuss the simplest model with a scale-invariant critical wave function, the quantum 
dimer model at the RK point. Here we also introduce the quantum Lifshitz model, the effective field 
theory of these new quantum critical points. In section II I II we generalize the relation between the 
quantum dimer model and the scalar field theory discussed in section ^ to include perturbations 
which drive the system in to a quantum disordered/deconfined phase or to a ordered/confined phase. 
In section IIVI we define the quantum eight-vertex model by finding a Hamiltonian whose ground-state 
wave function is related to the classical eight-vertex model. This will allow us to find quantum critical 
lines with variable critical exponents separating a Z2-ordered phase from a topologically-or dered phase . 
It will also allow us to place a number of previously- known models, in particular that of Ref. |43|, in 
a more general setting. We show in detail how to use the known results from the Baxter solution of 
the classical model to map out the critical behavior of the quantum theory. In particular we analyze 
in detail the confinement and deconfinement properties of the different phases and at criticality. In 
section O we study the non-abelian case, and see that the strongly-coupled limit of Yang- Mills theory 
with a Chern-Simons term has a wave functional local in two-dimensional classical fields Thi; 
theory is in a phase with topological order. In three appendices we give details of the correlators of 
the quantum Lifshitz field theory (Appendix^, and of the gauge-theory construction of the quantum 
six- vertex (Appendix IE)) and eight- vertex (Appendix [HI) models. 



II. SCALE-INVARIANT WAVE FUNCTIONS AND QUANTUM CRITICALITY 



The simplest lattice model we discuss is the quantum dimer model; the simplest field theory we dub 
the quantum Lifshitz theory. They both provide very nice illustrations of the properties discussed in 
the introduction. In the quantum dimer model, the space of states consists of close-packed hard-core 
dimers on a two-dimensional lattice. A quantum Hamiltonian therefore is an operator acting on this 
space of dimers, taking any dimer configuration to some linear combination of configurations. In every 
configuration exactly one dimer must touch every site, so any off-diagonal term in the Hamiltonian 
must necessarily move more than one dimer. The simplest such operator is called a "plaquette flip" : if 
one has two dimers on opposite sites of one plaquette, one can rotate the dimers around the plaquette 
without effecting any other dimers. For example, for the ith plaquette on the square lattice one has 




The operator Fi is defined as zero on any other dimer configuration around a plaquette (i.e. if the ith 
plaquette is not flippable). We define the operator Vi as the identity if the plaquette is flippable, and 
zero otherwise. 

The Rokhsar-Kivelson Hamiltonian for the quantum dimer model ^3] 

Hrk = Y^iVi - Fi) (2.2) 

i 

has the remarkable property that one can find its ground states exactly. They have energy zero, and 
every state (in a given sub-sector labeled by global conserved quantities) appears with equal amplitude 
in its ground-state wave function. These properties follow from the facts that HptK is self-adjoint, and 
(Vi — Fi)'^ = 2{Vi — Fi). Hamiltonians of the form H = QIQi necessarily have eigenvalues E obeying 
E > 0. Moreover, if one can find a state annihilated by all the Qi, then it is necessarily a ground 
state. The equal-amplitude sum over all states is indeed such a state. In the Schrodinger picture, the 
wave function for this state is easy to write down. Define Z as the number of all dimer configurations 
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in some finite volume. This is precisely the classical partition function of two-dimensional dimers with 
all configurations weighted equally. Then the properly-normalized ground-state wave function for any 
basis state \C) in the Hilbert space (i.e., any classical dimer configuration C) is 



(2.3) 



The wave function of the quantum system is indeed related to the classical system. 

One can extend this sort of analysis to compute equal-time correlators in the ground state. One finds 
simply that these correlators are given by the correlation functions of the two-dimensional classical 
theory. Thus for Hrk for dimers on the square lattice, one finds algebraic decay of the correlation 
functior is llOll . This model is then interpreted as a critical point between two ordered phases of the 
dimers |l3ll23i |. However, for the analogous Hamiltonian on the triangular lattice, the classical two- 
dimensional correlators are exponentially decaying |^ . One can also show that spinon-type excitations 
(sites without a dimer) are deconfined on the triangular lattice jlH. This means the quantum dimer 
model with Hrk on the triangular lattice is interpreted as being in a "liquid" phase, which has a mass 
gap and exponential decay of interactions, but which has no non-zero local order parameter. 

Such a relation between two-dimensional quantum theories is not limited to lattice models, nor are 
the ground-state wave functions required to be equal-amplitude sums over all configurations. We will 
construct now a simple (non-Lorentz invariant) two-dimensional quantum critical field theory i.e., a 
theory whose ground state wave function represents a two-dimensional conformal theory. 

Consider a free boson (/^(x, t) in two spatial dimensions and one time dimension. Instead of the 
usual Hamiltonian quadratic in derivatives, we use one which has been conjectured by Henley js^] to 
belong to the same universality class as the square-lattice quantum-dimer model. It is 



H 



where H = 99 as usual. The associated Euclidean action for the field (p is 

,2 



S 



(Fx 



(2.4) 



(2.5) 



This system, Eq. (|2.5() . also arises in three-dimensional classical statistical mechanics in the field-theory 
description of Lifshitz points 2], for example in (smectic) liquid crystals. For this reason we will call the 
system with Hamiltonian 1)2. 4() the quantum Lifshitz model Of particular relevance to our discussion 
is the long-ago observation by Grinstein [H^ ] that this system is analogous to the two-dimensional 
Euclidean free boson in that it represents a line of fixed points parametrized by k,. 

Let us rederive this result by quantizing the Hamiltonian 1)2. 4(1 . We impose the canonical commu- 
tation relations 



[<fix),U{x)]=i5{x-x') (2.6) 

so in the Schrodinger picture the canonical momentum is the functional derivative n(x) = —i6/6ip{x). 
The Schrodinger equation for the wave functional ^[f] is then 



d^x 



E^[ip]. 



Q(x) = i=(^ + .VV) , 



(2.7) 



We can find the ground-state wave function in the same fashion as we did for Hrk- Indeed, if we 
define 



(2. 
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the (normal-ordered) quantum Hamiltonian is then 

H = ^J cfx {q\x), Q(f)} - e,ac^ = J dPx Q\x)Q{x) (2.9) 

which is Hermitian and positive. Here V is the spatial volume (area) of the system, and we have 
normal-ordered the Hamiltonian by subtracting off the (UV divergent) zero-point energy density 

evac = - ^ lim Vl5{x - y) > . (2.10) 

Any state annihilated by Q{x) for all x must be a zero-energy ground state. The corresponding ground- 
state wave functional ([(/?] |vac) = "^oiv] satisfies Q^ofv'] = 0, where Q is defined in Eq. 1)2. 8|) . This is 
simply a first-order functional differential equation, and is easily solved, giving 



m^[^] = —e (2.11) 



where Z is the normalization 



-K / d^X {v^Y 
[Vif] e J . (2.12) 



The probability of finding the ground state in the configuration | [ip] ) is therefore 



-I -K d X {Vip) 
|^oMr = ^e J ■ (2.13) 

Consequently, the ground state expectation value of products of Hermitian local operators ©[(^(x)] 
reduces to expressions of the form 

I f —K / d'^x (V(/?)^ 

(vac|0[99(xi)] . . . 0[(^(x„)]|vac) = - / [V^] 0[^(xi)] . . . 0[(^(x„)] e J . (2.14) 

This two-dimensional quantum theory has a deep relation with a two-dimensional classical theory: 
the ground-state expectation value of all local observables are mapped one-to-one to correlators of a 
two-dimensional massless Euclidean free boson. The latter is a well-known conformal field theory, and 
its correlation functions are easily determined (for convenience we give them explicitly in Appendix 
This two-dimensional critical field theory is conformally-invariant, so the equal-time correlators of 
the quantum theory must refiect this. This scalar field theory is therefore not only quantum critical 
but it also has a time- independent conformal invariance. 

The equal-time expectation values of the "charge operators" 0[ip], as well as the correlation func- 
tions of the dual vortex (or "magnetic" ) operators discussed in Appendix exhibit a power-law 
behavior as a function of distance, as expected at a quantum critical point. As shown also in Ap- 
pendix their autocorrelation functions also exhibit scale invariance albeit with a dynamic critical 
exponent z = 2. This behavior of the equal-time correlator was shown earlier to occur in the quantum 
dimer model on the square lattice at the RK point: there is a massless "resonon" excitation, and the 
equal-time correlation functions for two static holons has a power law behavior equal to that of the 
monomer correlation function in the classical 2D dimer model on the square lattice. 

However, not all theories whose ground state can be found in this fashion need be critical with 
power-law behavior. As noted above, the quantum dimer model on the triangular lattice is not. As 
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we will discuss in section 11111 adding say a mass-like term to Q{x) in the scalar field theory gives a 
theory with exponentially-decaying correlations in the ground state. Thus one can understand phase 
transitions in such theories as well, and we will explore several of these in this paper. Nevertheless, we 
expect that when the quantum phase transition is continuous, the quantum critical points of generic 
theories of this type will have the basic structure of the quantum Lifshitz model. As discussed in the 
introduction, only quantum Lifshitz points can be conformal quantum critical points.^ 

In the next sections we will show that generalizations of this z = 2 quantum Lifshitz Hamiltonian 
also describe the quantum phase transitions between generalizations of the valence bond crystal states 
and quantum disordered states which describe deconfined topological fluid phases (provided these 
quantum phase transitions are continuous). Notice that phase transitions from deconfined to confined, 
uniform and translationally-invariant states are described by the standard Lorentz invariant z = 1 



These examples show that one can obtain precise information about some 2d quantum systems in 
terms of known properties of 2d classical systems. The trick of doing so is in finding a set of Qi (or 
Q{x) in the continuum) which annihilates the equal-amplitude state (or some other desired state), 
and then defining H = Y.iQ\Qi 0- This seems like it should be possible to do for any classical 2d 
theory, and indeed, there are many known examples of this sort. However, it is not clear for a given 
2d classical theory one can always find a Hamiltonian which is both local and ergodic (in this context, 
ergodic means that the Hamiltonian will eventually take the system through all of phase space with a 
given set of conserved quantum numbers). It is also not clear that even if such a Hamiltonian exists, 
whether it will have any physical relevance. 

Moreover, this simple relation of the ground-state wave function of a 2d quantum system to a 2d 
classical system is not at all generic: the quantum dimer model with Hhk and this quantum Lifshitz 
field theory are quite special. To illustrate this, let us discuss briefly the ground-state wave functions 
of standard quantum field theories at a (quantum) critical point. Consider first the most common case, 
the Lorentz invariant ip'^ field theory at criticality. Below L) = 4 space-time dimensions this critical 
theory is controlled by its non-trivial Wilson-Fisher fixed point. The resulting theory is massless and in 
general it has an anomalous dimension r/ 7^ 0. Scale and Lorentz invariance fully dictate the behavior 
of all the correlation functions at this fixed point. General fixed point theories are scale invariant and, 
in addition, they exhibit an enhanced, generally finite-dimensional, conformal symmetry. It is a very 
special feature of -D = 1 + 1-dimensional Lorentz-invariant fixed point theories that they exhibit a 
much larger, infinite-dimensional, conformal invariance. This enhanced symmetry leads to a plethora 
of critical behaviors in 1 -|- 1 dimensions. In contrast, there are relatively few known distinct critical 
points in higher dimensions for Lorentz-invariant field theories. 

It is well known that the knowledge of all the equal-time correlation functions determines completely 
the form of the ground state wave function, i. e. in the Schrodinger representation of the field theory |59|, 
IgoIIgiI ] . For a general theory, the ground-state wave function is a non-local and non-analytic functional 
of the field configuration. Thus, at the Wilson-Fisher fixed point, which describes theories with 
only a global conformal (scale) invariance, the structure of the ground state wave function is quite 
complicated. For instance, the probability of a constant field configuration (/^(x) = (/? at a critical point 
has the universal form |6l|| 



For a Lorentz-invariant Lp"^ theory the universal critical exponent is given hy 6 = {d + 2 — 'rf)/{d — 2 + r]), 
and A and B are two non- universal constants. Therefore, at criticality the wave function in general is a 
non-analytic non-local functional of the field configuration. In contrast, the ground-state wave function 




1+5 



(2.15) 



^ The role of conformal invari ance in 2D classical dynamics with z — 2 and anisotropic 3D classical Lifshitz points was 
considered recently in Refs. |55l56| . 
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of a 1 + 1-dimensional relativistic interacting fermions (a Luttinger liquid), which is a conformal field 
theory, has a universal non-local non-analytic Jastrow-like power law factorized form 'g^I consistent 
with the form found by the Bethe-ansatz solution of the Calogero-Sutherland model ^63..,64.65il] . This 
structure is a consequence of the (local) conformal invariance of the 1 + 1-dimensional theory. Even 
in 1 + 1 dimensions, the wave function is generically non-local. 



III. DIMERS, FERMIONS AND THE QUANTUM LIFSHITZ FIELD THEORY 

In Section ^ we showed how to find the exact ground-state wave functions of the quantum dimer 
model with Rokhsar-Kivelson Hamiltonian Hfix, and the quantum Lifshitz scalar field theory. In this 
section, we will describe their properties in more detail, and some simple generalizations. 



A. Prom the square to the triangular lattice 



The ground states of Hrk are the sum over all classical dimer configurations in a sector with equal 
amplitudes 1^. There are a number of useful generalizations to models where the ground state is still 
a sum over all states in a sector, but not necessarily with equal amplitudes. One interesting case is a 
quantum dimer model which interpolates between the square and triangular lattices. 

A triangular lattice can be made from a square lattice by adding bonds across all the diagonals in 
one direction. The classical dimer model on the square lattice can be deformed continuously into the 
triangular-lattice model by assigning a variable weight w for allowing dimers along these diagonals: 
for w = we have the original square-lattice model, while w = 1 gives the triangular-lattice one |ll|. 
There is a two-dimensional quantum Hamiltonian which has the w-dependent classical dimer model 
as its ground state. For dimers on opposite sides of a plaquette of the original square lattice, the 
Hamiltonian remains H^k- In addition, however, parallel dimers on adjacent diagonals can also be 
flipped see figure ^ Like the flip on the square lattice, this flip can be done without violating the 
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FIG. 1: Dimer flips on the triangular lattice 



close-packing and hard-core constraints. 

When the off-diagonal terms in the Hamiltonian consist solely of these flips, H^, breaks up into 
two-by-two blocks like Hjik does. In the ground state dimers along diagonals should get a weight w, 
so the Hamiltonian must explicitly depend on w. To construct a Hamiltonian with the desired 
ground state, we find a set of operators Qi, each of which annihilates this state. For example, denote 
a configuration with dimers on adjacent diagonals as |1), and |2) as the configuration to which it is 
flipped, as shown in fig. ^ The Qi acting on these two configurations is then 

Q^ = ^A—. (3-1) 
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where the first row and column correspond to state |1), while the second corresponds to state |2). 
There are three types of Qi: q^^*!^^^^) g^g^^g dimers on opposite sides of a plaquette of the original 
square lattice, while (^j^^^^^) and Q^^^^^^ are associated to flips involving dimers on the 'diagonal' 
links; see figure ^ The operator q^^^^^'^^) given by H3.1|) with w = 1. We then take 

H^ = Y. [qP^^^''^ + ^^r''^ + (3.2) 

i 

where the sum is over all plaquettes i. Each Qi is a projection operator, so Qi = Q\ = [Qi)'^ ■ 

We have defined the operator Qi so that it annihilates the state nP'\l) + |2). The ground state of 
Hyj is then the sum over classical dimer model states with each state weighted by , where T) is the 
number of dimers along diagonals in that state. More precisely, one must find the conserved quantities 
for a given value of w and boundary conditions; the sum over all states with weight nP in that sector 
is an eigenstate of Hrk with zero energy. The ground-state wave function for a configuration with T) 
diagonal dimers is then 

MV] = (3.3) 

where Z{w) is the partition function for the classical dimer model with diagonal dimers receiving weight 
w. Note the nP' in the argument of Z in the denominator: this is because in quantum mechanics 
probabilities are given by l^'oP- For the square or triangular lattice, this is unimportant, because 
u; = or It; = 1, both of which have = w. For = 1, we have the equal-amplitude sum over 
all dimer states of the triangular lattice, the model discussed in Ref. [^. For u) = 0, we recover a 
slight generalization of the original square-lattice quantum-dimer model of Rokhsar and Kivelson ^3] ■ 
In this limit, this Hamiltonian reduces to H^ik plus a potential term forbidding dimers on adjacent 
diagonal links. Isolated diagonal dimers are still allowed for w — > 0, but since none of them can be 
fiipped, the Hamiltonian does not affect them at all. Thus they can be viewed as fixed zero-energy 
defects in the square-lattice quantum-dimer model. The ground-state wave function for a given set of 
defects is the equal-amplitude over all configurations of dimers on the sites without defects. 

Since we know the exact ground-state wave function for any one would like to compute the 
correlation functions in this model. In most two-dimensional lattice models, even those solvable by 
the Bethe ansatz, this is extraordinarily difficult or impossible. However, the classical dimer model 
is special in that one can do such computations, because like the two-dimensional Ising model, it is 
essentially free-fermionic. Precisely, its partition function and correlators can be written in terms of 
the Pfaffian (the square root of the determinant) of known matrices [g^. One can rewrite the Pfaffians 
in terms of a functional integral over Grassmann variables at every site on the lattice [g^I ■ The action 
in the case of equal Boltzmann weights is quadratic in the Grassmann variables, so one can compute 
easily any ground-state correlation function using the dimers, because the dimers can be written in 
terms of the fermions. This was discussed for the triangular lattice in ;8]. The correlators of spinon-like 
or holon-like excitations are much more complicated, but the computation was done for w = in j6Q|] , 
and for arbitrary w in On the lattice, the holon is a defect or monomer, a site without a dimer. 
The holon-creation operators are not local in terms of the fermionic variables, in a manner reminiscent 
of how the spin and fermion operators are non-local with respect to each other in the two-dimensional 
Ising model. The holon two-point function is valuable in that it gives an order parameter for the 
phase with topological order: if it is non-vanishing as two holons are taken far apart, the holons are 
deconfined and we are indeed in a topological phase. The existence of topological order was previously 
established for the triangular lattice w = in 1^3] the explicit correlator was computed, and 

indicates the topologically-ordered phase exists for any non-zero w. 

The ground state of the quantum dimer model with Hy^ is therefore well understood for any w. 
There are no exact results for the excited states, however. In fact, since Hyj does not have any action 
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upon empty sites, one can give the holon any gap desired without changing the ground states. It is 
therefore useful to find a continuum hmit and study the field theory describing this model. In other 
words, we would like to understand a field theory with partition function equal to the continuum 
limit of partition function Z(w) in Eq. (|3.3|) . Since ground-state correlators for the square lattice are 
algebraically decaying, the w = model is critical and should have a sensible continuum limit. Indeed, 
when w = 0, the Grassmann variables turn into a single free massless Dirac fermion field; its action is 
the usual rotationally-invariant kinetic term. The dimer correlation length ^ was computed exactly as 
a function of w; for the triangular lattice it is about one lattice spacing, while ^ diverges as 1/w for w 
small |llt|. Thus there is a field-theory description of the continuum limit of the classical dimer model 
valid as long as w is scaled to zero with the lattice spacing a such that w/a remains finite. Since the 
action is still quadratic in the Grassmann variables, the resultin g fe rmionic field theory remains free. 
However, the Dirac fermion receives a mass proportional to w/a 11]. 

There are several useful aspects of taking the continuum limit, apart from finding the excited-state 
spectrum. Gorrelators are easier to compute: for w = one can use conformal field theory [i^, while 
in the scaling limit u; — > one can use form-factor techniques |^ . Another useful fact is that (ignoring 
boundary conditions), a Dirac fermion can be described in terms of two decoupled Ising field theories. 
In the continuum, the holon can be written in terms of the product of the spin field in one Ising model 
with the disorder field in the other Ising model. Taking w away from zero amounts to g iving one Ising 
order field an expectation value, and the other disorder field an expectation value. [ll|. Thus one can 
see directly in continuum that the holon order parameter is non- vanishing for w ^ 0. We will see in 
the next section that one can also understand the physics of the quantum Lifshitz critical line for all 
K in terms of two (coupled) Ising models. 



B. The critical field theory 

We would therefore like to find a natural-looking quantum field-theory Hamiltonian which has as 
its ground-state wave functional 

g — SoiracH'] 

V ^ Dirac 

where Soirac is the usual action for a rotationally-invariant action for a free Dirac fermion in two 
Euclidean dimensions. Since this wave functional involves Grassman numbers, the easiest way to 
think of l^oP as a weight in the path integral defining all correlators. While it is possible to find 
a Hamiltonian acting on this fermionic basis, it is more convenient and more intuitive to instead 
use bosonic variables. It is more convenient because correlators in a massless Dirac fermion theory 
(including those involving the product of spin fields) can be bosonized, meaning that they can be 



written in terms of correlators of free scalar fields |48ll69l | . It is more intuitive because the classical 
dimer model on the square lattice has a simple description in terms of a "height" variable [s^]. A 
height is an integer-valued variable, which typically in the continuum limit turns into a scalar field. 
This description will allow us also to make contact with the quantum Lifshitz model. 

Recently, Moessner et al. |0] generalized Henley's argument of Ref . [s^] , and used the connection 
between quantum dimer models and their dual quantum roug hening (height) models [l3|l3|23|730, 
It^It^It^ I to argue that the Hamiltonian of Eq. ()2.4|) actually defines the universality class of quantum 
critical points between valence bond crystal phases. The nature of the phase transition between valence 
bond crystal states is the focus of much current research. Quite recent results by Vishwanath et al. [t^ 
^, and by Fradkin et al. 76:], show that the transition between valence bond crystals is generically first 
order, as expected from a simple Landau argument. Nevertheless, when the transition is continuous. 



^ We became aware of Ref. as this paper was being completed. 
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it is described by the Hamiltonian of Eq. (|2.4j) which must be regarded as a (rather rich) multicritical 
point. 

To see how the height description arises, let us go back to the square-lattice quantum dimer model, 
where reduces to Hrk- To map the square-lattice classical dimer model onto a height model, 
one first assigns a height variable to each plaquette. In going around a vertex on the even sub-lattice 
clockwise, the height changes by -|-3 if a dimer is present on the link between the plaquettes, and by 
— 1 if no dimer is present on that link. On the odd sub-lattice, the heights change by —3 and -|-1 
respectively. The flip operator Fi on a plaquette i changes the height on that plaquette by either 
±4. To take the continuum limit, is convenient to turn this into a model with heights on the sites. 
We define h on each site to be the average value of the four plaquette heights around that site^; see 
figure 121 To avoid overcounting configurations, we identify the height h with h + ^. The flip operator 
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F corresponds to changing /i — > /i it 1 on all four sites around a plaquette (the it depending on the 
sub-lattice). Columnar order for dimers corresponds to an expectation value for /i, while staggered 
order for dimers corresponds to an expectation value for dh. One can obtain a Hamiltonian with 
ordered ground states by allowing the coefficients of the two terms in Eq. 1)2. 2() to be different. If the 
coefficient of Fi is larger, this favors columnar order; if the coefficient of is larger, staggered order is 
favored. It is widely assumed that Hr^^, where the coefficients are equal, describes a phase transition 
between the two kinds of order 12, 23,70 "t^. However, this has never been proven, and there exists 
the possibility of an intermediate "plaquette" phase F^. 

We would now like to take the continuum limit of the quantum dimer model in its height description. 
In this limit, we identify the height h with a scalar field 4:{p{x). Like the height, the scalar field must 
be periodic, so we identify (p with (p + 1. We have already noted that the continuum correlators of 
the square-lattice quantum dimer model are those of a massless Dirac fermion. The correlators in 
the ground state of the quantum Hamiltonian in terms of <p must be identical. The quantum Lifshitz 
Hamiltonian Eq. ()2.4() and Eq. (|2.9|) has ground-state correlators of the form given in Eq. 1)2. 14() . When 
= 27r, these correlators are precisely those of a Dirac fermion; this is a result of the widely-known 
procedure known as bosonization ji^. In fact, correlators of many two-dimensional critical classical 
statistical mechanical systems, not just free fermions, can be written in terms of exponentials of a 
free boson |6^; we collect some of these results in Appendix In the next section we will display 
quantum lattice models whose continuum limit corresponds to all values of k. 

The Hamiltonian Hhk on the square lattice in the continuum limit is therefore identified with 
the quantum Lifshitz Hamiltonian with = 2tt. This also allows a qualitative understanding of 
the physics away from the RK point [s^]. A phase with staggered order should have an expectation 
value of dx^ or dyip in the continuum theory. Adding a term (Vip)'^ to the Hamiltonian with negative 
coefficient will drive the system into such an ordered phase. Adding this term with positive coefficient 
will favor a constant value of ip, driving the system into columnar order. Adding terms like cos(27r(/3) 
will also drive the system into a phase with columnar order. Thus one expects that at a critical 
point like that described by Hrk on the square lattice, the coefficients of (Vip)'^ and cos(2n7r</j) will 



* This construction has been employed extensively in statistical mechanics models of classical dimers and loops; see for 
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vanish. This leaves Eq. (|2.4j) as the simplest non-trivial Hamiltonian with the desired properties. The 
requirement that the ground state be equivalent to a free fermion then fixes the coefficient k; note 
that if desired n can scaled out of the Hamiltonian by redefining the compactification relation to be 

(f (f + y/K. 

It is not at all clear whether critical 2 + 1-dimensional field theory and the continuum limit of 
Hrk are identical for excited states, although the above heuristic argument is very suggestive. The 
Hamiltonian for the scalar field theory Eq. 1)2. 9() is purely quadratic in the field (p, so one can obtain 
essentially any desired information exactly. Hjik is not so simple on the lattice, but one may hope the 
two models are in the same universality class. In any event, we can easily extract all the excited-state 
energies for the quantum Lifshitz field theory. The operators Q{x) and Q^{x) are essentially harmonic- 
oscillator creation and annihilation operators: the equal-time commutation relation Eq. ()2.6|) implies 
that 

Q(x),Qt(y)] =KVV2)(x-y) (3.5) 

The ground state is indeed annihilated by all Q{x), so we can create excited states by acting with 
Q^{x). The commutation relations Eq. ()3.5|) mean that the dispersion relation is E = Kp^. This theory 
is gapless but not Lorentz- invariant: the dynamical critical exponent is z = 2. 

The exponent z = 2 can also be seen by looking at the classical action associated with this Hamilto- 
nian. This will also allow us to make contact with the earlier three-dimensional statistical-mechanical 
results. The action consistent with the Hamiltonian of Eq. 1)2. 4|1 and the canonical commutation 
relations of Eq. 1)2. 6|) is 



d^x 



(3.6) 



Clearly, this action is not Lorentz invariant and has z = 2. It is rotationally invariant only in the XY 
plane. Defining the imaginary time t = it gives the Euclidean action ()2.5() . The imaginary time axis 
r can be regarded as the z-coordinate of a three-dimensional classical system in which (p{x, r) is an 
angle-like variable and the action represents the spin-wave approximation of an anisotropic classical 
XY model. In general one would have expected a term proportional to the operator (Vip)'^ with a 
finite positive stiffness in the plane. This is so in the XY ferromagnetic phase. On the other hand, if 
the stiffness becomes negative, there is an instability to a modulated helical phase. The action of Eq. 
()2.5() represents the Lifshitz point, the critical point of this phase transition 54,8o|] where the stiffness 
vanishes. This effective action also plays a central role in the smectic A-C transition 81] and in other 
classical liquid crystal phase transitions associated with the spontaneous partial breaking of translation 
and/or rotational invariance (2ii82l] . The compactified version of the problem {i.e., the identification 
(f {p -\- 1) has also been considered in this context, for example in Ref. 54]. The choice of period in 
general depends on the physical context of the problem. 

The square-lattice quantum dimer model and the scalar field theory with Hamiltonian Eq. (|2.9|) 
are both at critical points, in that the correlators in the ground state are algebraically decaying. The 
quantum dynamics implied by this Hamiltonian must be compatible with the 2D time-independent 
conformal invariance. In particular, the spectrum of the quantum theory must be gapless and, as we 
learned from this example, the dynamic critical exponent must he z = 2. Notice however, that z = 2 
alone does not guarantee a gapless (or even critical) theory. Indeed, instructive counter-examples to 
this statement are well known in the theory of (the absence of) quantum roughening [s^ls^ where 
quantum fluctuations destroy the critical behavior and lead to an ordered state through an order-from- 
disorder mechanism. 
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C. The off-critical field theory 

We have thus shown that the continuum hmit of the square-lattice quantum dimer model is de- 
scribed by the quantum Lifshitz model at a special point = 2tt. We also argued that at least 
some deformations of the two models result in ordered phases. However, we saw at the beginning of 
this section that not all deformations of the square-lattice quantum dimer model result in an ordered 
phase. Allowing dimers across the diagonals with Hamiltonian Hyij results in a topologically-ordered 
phase, where the order parameter is not local. 

In this subsection, we find a bosonic field theory describing the topological phase in the continuum 
limit. We showed above that in the scaling limit — > with w/a finite, the ground-state wave function 
Eq. (|3.4[) can be written in terms of a free massive Dirac fermion of mass proportional to w/a. In two 
dimensions, the bosonic version of a massive Dirac fermion is the sine-Gordon model at a particular 
coupling. Precisely, the two-dimensional fermion action is equivalent to 

= J (fx [K{Viff - Acos(27r99)] . (3.7) 

For the free- fermion case, we have n = l/27r; we will discuss the more general case in section ITVl In the 
fermion language, different values of k correspond to adding a four-fermion coupling to 5*2^. To find a 
Hamiltonian with this two-dimensional action describing the ground state, we again find an operator 
Q{x) and define the Hamiltonian via Eq. ()2.9() . The operator 

Q(x) = i= (^^ + KVV + ^sin(2^(^)^ (3.8) 

annihilates the wave functional \I' cx e~^'^'^. Because of the extra term in Q, the commutator 
[Q{x), Q^{y)] is not a simple c-number, but in fact depends on the field configuration 99, 

Q{x),Q\y)\ = kV^ (5(2)(x-y)) + Asin(27rvp(x))5(2)(f -y-) . (3.9) 

Thus, normal-ordering the Hamiltonian is not just an innocent ground state energy shift: the two 
parts of (|2.9|) are not the same here. To obtain the desired ground-state wave functional, we must 
define H of the form f Q^Q. The three-dimensional version of this model was discussed in Ref. [s^ . 

This Hamiltonian is not quadratic in the field (p except at the critical point m = 0, so that even with 
this fine tuning this model cannot be solved simply. Since Q and do not have simple commutation 
relations, we cannot simply find the spectrum of this theory. Of course, one can compute properties 
in the fermionic picture, but as noted before, computations involving spin fields are non-trivial in this 
basis as well. However, since there are dimensionful parameters in the Hamiltonian and no spontaneous 
breaking of a continuous symmetry, it seems likely that the Hamiltonian is gapped. Moreover, in the 
limit with X/k finite, the action Eq. 1)3. 7jl reduces to that of a free massive boson. Then one can solve 
the model explicitly, and the quantum Hamiltonian indeed has a gap. When A is reduced to a finite 
value, the gap should remain^. 



^ The equal-time correlators of the vertex operators defined and computed in Appendix EI can also be computed away 
from the critical point by either by a naive semi-classical argument, which predicts the simple exponential decay we 
just discussed, or in a more sophisticated way by means of form factors [g^I . Apart from a subtle bound state structure 
in the spectral functions, the more sophisticated approach confirms the essence of the naive semi-classical result. On 
the other hand, the z = 2 character of the critical theory suggests that time-dependent correlation functions must be 
consistent with this fact, and that the correlators must be functions of or t, and that the time dependent Euclidean 
auto-correlation functions may obtained from the equal-time correlator by replacing x'^ ^ \t\. These arguments are 
consistent with the renormalization-group results of Ref. [s^ 
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There are many terms in the Hamiltonian of this field theory, and their coefficients must be fine- 
tuned to enable us to compute the ground-state wave function explicitly. There are terms like cos(47r(/j) 
and (Vip)^ which, as noted above, tend to order the system. However, the exact lattice results for 
from |ll| show that the ground-state correlators for w small are those of the bosonic Hamiltonian with 
these special couplings. Thus this model is not ordered but rather topologically ordered. Moreover, 
since the model is gapped, we expect that its physical properties are robust and persist even when 
the coefficients are tuned away from this special point. Thus there must exist a topological phase, not 
just an isolated point. An interesting open problem is to understand how large (in coupling constant 
space) the topological phase is, as compared to the ordered phases. 



IV. A QUANTUM EIGHT- VERTEX MODEL 

In the last section, we discussed a 2 + 1-dimensional theory whose ground-state wave function is 
simply described in terms of a classical two-dimensional bosonic field. With vanishing potential and a 
particular value of the coupling k, it is believed to describe the continuum limit of the quantum dimer 
model on the square lattice. In this section, we will study lattice models which in the continuum limit 
allow arbitrary values of k. These models also have a quantum critical line separating an ordered 
phase from a topologically-ordered phase. 

The degrees of freedom in our model are those of the classical two-dimensional eight-vertex model. 
These are arrows placed on the links of a square lattice, with the restriction that the number of arrows 
pointing in at each vertex is even. This means that there are eight possible configurations at each 
vertex, which we display in figure (jSJ. The classical Boltzmann weights for a given vertex in the 




a a b be c d d 

FIG. 3: The eight vertices and their Boltzmann weights 



zero-field eight-vertex model are usually denoted by a, b, c and d, as shown in the figure. Since we 
are interested in rotationally-invariant theories, we set a = 6 in the following; moreover, since we can 
rescale all the weights by a constant, we set a = b = 1. A typical configuration is displayed in figure 
(jlj ; the Boltzmann weight of such a configuration is given by the product of Boltzmann weights of the 
vertices. 

The classical eight-vertex model is integrable, and many of its properties can be derived exactly 
(SS*]. For a = 6 = 1, it has ordered phases for c > d + 2 and d > c + 2. In these phases the Z2 symmetry 
of flipping all the arrows is spontaneously broken. Critical lines with continuously varying exponents 
at c = d + 2 and d = c + 2 separate the ordered phases from the disordered one |c — (i| < 2. The 
correlation length diverges as |85| 

e~ ||c-d| -2p^/^^''\ ;u = 2tan-i(y^) , (4.1) 

near these critical lines. (For vr/^u an even integer this is multiplied by log||c — d\ — 2\ .) When 
c = 0, d < 2 or d = 0, c < 2, the model is also critical; in fact the partition function on this line 
can be mapped onto that for the order-disorder critical line. For d = 0, the exponent in Eq. (|4.1|) 
diverges: there is a Kosterlitz-Thouless transition as one brings c through 2. Another useful result for 
the classical correlation length is that it is zero on the line c = d; this is the state of maximal disorder. 
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FIG. 4: A 



iypical configuration in the eight-vertex model 



An order parameter which will be useful later comes by rewriting the model in terms of an Ising 
spin at the center of each plaquette. This description is best thought of as two Ising models, with spins 
t{A) on one sublattice, and t{B) on the other. Then the Boltzmann weights can be written in terms of 
two Ising couplings between nearest sites on the A lattice and on the B lattice, and a four-spin coupling 
between the two A and two B spins around a site of the original lattice. The polarization operator 
of the eight-vertex model becomes t{A)t{B). One finds that its expectation value is non-vanishing in 
the ordered phase, and vanishes in the disordered phase |c — d| < 2 [s^. One also can define Neel-like 
staggered order parameters, in terms of t{A) and t{B) individually, which do not vanish in the ordered 
phase. Along the line cd = 1, the four-spin coupling vanishes, so the eight-vertex model turns into two 
decoupled Ising models; the model here in this case be solved by using Pfaffian techniques js^. Along 
the line c = d, the two Ising couplings vanish, leaving only the four-spin coupling. Thus this line in 
the classical model has an extra Z2 gauge symmetry. 

The classical eight- vertex model has a number of useful dualities '85]. They can be described by 
defining the combinations Wi = {a + b)/2, W2 = {a - b)/2, W3 = (c + d)/2 and VF4 = (c - d)/2. The 
partition function is invariant under the exchange of any two of the Wj and under the Wj — Wj for 
any j. These dualities, for example, map the critical line c = d + 2 to the critical line d = 0, c < 2 by 
exchanging Wi with W4^. In Ising language this amounts to performing Kramers- Wannier duality on 
one of the two types of Ising spins. The line c = d + 2 is invariant under the exchange Wi VF3. In 
Ising language, this duality amounts to taking the Kramers- Wannier dual of both types of Ising spins. 
Denoting the dual spins as n{A) and fJ,{B), duality means therefore that in the disordered phases, the 
expectation value {t{A)t{B)) is non-vanishing. 



A. Construction of the quantum eight-vertex Hamiltonian 

We now define a quantum Hamiltonian acting on a Hilbert space whose basis elements are the 
states of this classical eight-vertex model. To define such a Hamiltonian, we first need the analog of 
the flip operator in the quantum dimer model. A flip operator needs to be ergodic: by flips on various 
plaquettes one should be able to reach all the states with the same global conserved quantities. The 
simplest such operator for the eight-vertex model is the operator which reverses all the arrows around 
a given plaquette. We write this flip operator J^i explicitly in gauge-theory language in appendix lO 
Note that as opposed to the quantum dimer model on the square or triangular lattice, all configurations 
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in the quantum eight- vertex model are flippable : J-i preserves the restriction that an even number of 
arrows be pointing in or out at each vertex. 

The simplest Hamiltonian has no potential energy, just a flip term. It is convenient to write this 
in terms of a projection operator: for / the identity matrix, we have (/ — !Fi)'^ = 2{I — Ti). Then the 
Hamiltonian 



has a ground state corresponding to the equal-amplitude sum over all eight-vertex model states. In 
terms of the Boltzmann weights introduced above, this is the state with a = h = c = d = \. A 
Hamiltonian with the same ground state was introduced by Kitaev There the eight-vertex- model 
restriction of having an even number of arrows in and out at each vertex was not required a priori, but 
instead a term was introduced giving a positive energy to vertices not obeying the restriction. The 
zero-energy ground state can therefore include no such vertices, so the ground state for the model of 
is indeed the sum over the eight- vertex-mo del configurations with equal weights. Because every 
plaquette is flippable and all configurations have equal weights, different Ti commute. Therefore all 
the terms in Eq. 1)4. 2() commute with each other, so they can be simultaneously diagonalized and their 
eigenstates can easily be found, as for the quantum Lifshitz field theory. The model is gapped, and is 
in a topologically-ordered phase ;42] . This follows as well from the results for the classical eight-vertex 
model discussed above: for c = d = 1 all Ising couplings vanish resulting in two decoupled Ising models 
at infinite temperature. At this point the order parameter vanishes, but the non-local order parameter 
does not. We can thus interpret the product of dual Ising variables fi{A)fj,{B) as a topological order 
parameter. 

This ground state can also be mapped onto the ground state of a Z2 gauge theory deep in its 
deconfined phase. In Appendix O we give a detailed derivation of the Z2 gauge theory of the full 
quantum eight-vertex model we define below, as well as the its dual theory which we will use to 
characterize some of the phases. In Appendix^ we discuss the U{1) gauge theory description of the 
quantum six-vertex model, which describes the limit d = 0, and by duality, the lines c = 0, = + 2 
and d"^ = + 2. 

We now find a two-parameter quantum Hamiltonian whose ground state is a sum over the states of 
the eight-vertex model with amplitudes given by the classical Boltzmann weights with arbitrary c and 
d. We are still keeping a = 6 = 1 to preserve two-dimensional rotational invariance, but this restriction 
can be relaxed if desired. Our model is neither the simplest nor the most natural extension of the 
classical eight-vertex model: a simpler Hamiltonian was proposed by Chakravarty [s^ in the context 
of d-density waves. Although we believe that our Hamiltonian and Chakravarty's describe the same 
physics, we are not aware of any simple mapping between these models. Another lattice model related 
to the quantum eight- vertex model discussed here was introduced in Ref. 'QQl. The main virtue of the 
construction that we use here is the structure of the ground-state wave function. 

Finding a quantum Hamiltonian with a known ground state is straightforward to do by using the 
trick discussed above (and in Ref. [HH). Namely, we find a Hamiltonian of the form 



where Qj = Q| tx Q?. To yield the desired ground state, each operator Qi must annihilate the sum 
over states with each state weighted by c^'^d^'^, where Nc and are the number of c and d type 



However, the quantum dimer models on the Kagome j87| and Fisher |88| lattices have Hamiltonians where all plaquettes 






(4.3) 



are flippable. 
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vertices in that state. In particular, we look for a Qi of the form 



2^ = E[ 



(4.4) 



where Vi is diagonal and depends on the Boltzmann weights for the four vertices at the corners of the 
plaquette i. Since {J-i)^ = /, this Hamiltonian breaks into 2 by 2 blocks like Hrk for the quantum 
dimer model. If we choose the potential Vi so that the blocks are of the form 



V -1 
-1 



-1 ) , (4.5) 



Qi will have the desired properties. 

To find Vi, let ric be the number of c-vertices at the corners of the plaquette i, and let ric be the 
number of c-vertices around the plaquette after it is flipped by J-i. Likewise, let rid be the number 
of d-vertices around the plaquette, while Ud is the number of d vertices in the flipped configuration. 
Note that Ti always flips a c or d vertex to an a or 6 vertex, and vice versa. Consequently, we take 
the operators Vi to be of the form 

Explicitly, we can write the projectors as follows 

\ —1 (fic-n^ ffii-n,i j ■ {^■') 

This Qi is indeed proportional to a projection operator: the Hamiltonian (|4.3() has the classical eight- 
vertex model as an eigenstate. This holds for any choice of the Wi, but in the Z2 gauge-theory language 
of appendix[nj it is natural to set all Wi = 1. The amplitude of the ground-state wave function of this 
HqSv on a state with Nc c-vertices and d-vertices is 

where Z(c, d) is the partition function of the classical two-dimensional eight- vertex model with weights 
a = 1,6 = 1, c and d. The arguments in the denominator are (? and because averages in the quantum 
model are calculated with respect to l^'oP- 

In section lllll we discussed how by taking the limit tt; ^ in i?^, one can recover a slightly- 
generalized square-lattice quantum dimer model which allows defects with no dynamics. Similarly, 
here one can find a quantum six-vertex model with some defects allowed by taking the limit d — > 
0.'' Configurations with > rid on any plaquette will receive infinite potential energy and so are 
disallowed. Some useful facts are that nc + nd + Uc + nd = 4, and ric — ric = Q mod 2 and rid — rid = 
mod 2. When d — > 0, configurations with = = are the flippable plaquettes of the six- vertex 
model, and are obviously included in the ground-state wave function. Configurations with = = 1 
or rid = rid = 2 are not suppressed, but the flip therefore preserves the number of d vertices around each 
plaquette. Thus in the d — > limit we can view these d vertices as defects in the quantum six-vertex 
model. The ground-state wave function is a sum over all allowed states with a given set of plaquettes 
with defects. The amplitude of each configuration in this sum is proportional to c^''. We discuss the 
relation between the quantum six- and eight-vertex models in more detail in the Appendices IHl and lUl 



^ This Hamiltonian in the special case c = 1, d = was discussed in Ref. Moreover, the quantum six- vertex model 
should be in the same universality class as the"supersymmetric" XY model introduced some time ago ;92i; this model 
defines a similar Hamiltonian acting on the 2D classical XY model, and presumably can also be understood in terms 
of a mapping to the quantum Lifshitz model. Another quantum six-vertex model has been proposed as a model of a 
planar pyrochlore lattice in Ref. [o^ . This six- vertex model is not of the Rokhsar-Kivelson type; it is the six- vertex 
limit of the simpler quantum eight-vertex model we discuss in section HV Dl 



19 



Ordered 



Dual 6- vertex 




Dual 6-vertex 



6- vertex 



2 



,2 



FIG. 5: Phase diagram of the quantum eight-vertex model: phases / and // are separated by a dual 6-vertex 
transition (same with I and ///); 6-vertex denotes the 6-vertex model critical lines and KT are 2D Kosterlitz- 
Thouless transitions; the dotted line shows that the Kitaev point is smoothly connected to the critical regime 
of the eight-vertex model. 



Now that we have the exact ground-state wave function for the quantum eight-vertex Hamiltonian 
of ()4.3(1 . we can use it to determine the phase diagram. Given the fact that the probability of a 
configuration of arrows in the ground state wave function is equal to the Boltzmann weight of a 
classical eight-vertex model, we can deduce much of the physics of the quantum theory (at least its 
equal-time properties) directly from the Baxter solution of the classical 2D eight-vertex model js^, as 
well as from Kadanoff's classic work on its critical behavior [6^l93l95j |. The only change is that the 
weights must be squared here, since in quantum mechanics we weigh configurations with The 
phase diagram is displayed in figure El note that the axes are labeled by and cP. 

We will now use this knowledge, as well a simple perturbative arguments in the quantum theory, to 
determine the phase diagram, the behavior of physical observables in the different phases, and (much 
of) their critical behavior. A useful fact is that for cd = 1, the classical model with partition function 
Z{cP, 1/c^) is equivalent to two decoupled Ising models. This decoupling is a property only of the wave 
function, not of the Hamiltonian of the full 2 -|- 1-dimensional quantum theory. Since this line cd = 1 
goes through both ordered and disordered phases, much of the physics of the quantum eight-vertex 
model can be described at least qualitatively in terms of decoupled Ising models. In particular, for any 
values of c and d except those on the critical line, correlators decay exponentially fast with distance 
with a correlation length which diverges as the phase boundary is approached, in a manner given in 
Eq. ()4.1|) . This exponential decay occurs in general, not just on the decoupling curve. 

As with the dimer models discussed in section HTll the partition function for cd = 1 can be expressed 
in terms of Grassmann variables with only quadratic terms, i.e. free fermions. Duality means that 



B. The phase diagram 
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FIG. 6: The ordered phase "antiferroelectric" or "staggered flux" phase of the eight- vertex model for c ^ 1. In 
this limit there are only c vertices in the ground state. 



Z^c^jd?) is free-fermionic for + d'^ = 2 as well. The correlators in the continuum limit of the 
(critical) square-lattice quantum dimer model are therefore identical to those obtained for c = \/2, 
c? = 0. The special point c = d = 1 discussed above and in Ref. 143] (labeled in fig. [51 as "Kitaev") 
is also free fermionic. At this point one does not need the Pfaffian techniques to compute correlators 
exactly, and one finds that the model is in a disordered phase in the Ising-spin language. However, 
the expectation value {^{A)^{B)) is non- vanishing, so there is topological order at this point. 
Let us now discuss the different phases of this system. 

1. The Ordered (Confined) Phase: 

From the known phase diagram of the classical eight-vertex model [s^, we conclude that the 
ground state of the quantum model with the Hamiltonian of Eq. (|4.3|1 has an ordered phase 
for (? > + 2 (and also for d? > + 2). That this is an ordered phase can be seen easily 
by considering the limit c ^ oo (with d fixed). In this limit the ground state is dominated 
by just two configurations, related to each other by a lattice translation of one lattice spacing, 
which have a c vertex on every site, as shown in fig. [HI In this phase the staggered polarization 
operator {t{A)t{B)) has a non-vanishing expectation value. This result can also be obtained 
directly from the Hamiltonian of the Z2 gauge theory, Eq. (|4.1I|) . since for c large the potential 
energy term Hy dominates and in it the piece associated with the c projection operators. 

The ordered phase is also confining. Below we will discuss the behavior of the Wilson loop 
operator and show that in this regime it obeys an area law, the hallmark of confinement |3l]l57| |. 
We will also show that the energy of a state with two static sources grows linearly with their 
separation. We should also note here that the equal-time fermion correlation function has an 
exponential decay in this phase, suggesting that this phase may support massive fermionic 
excitations. 

2. The Disordered (Deconfined) Phase: 

From the exact solution of the classical model, we know that there is a disordered phase for 
< d'^ + 2 (d^ < + 2) (see Fig. This is also a deconfined phase. This is most easily 



seen by taking a point deep in this phase, such as the Kitaev point a = h = c = d = l [4j]. In 
fact, all along the line a = 6 = 1 and c = d, the correlation length of the eight-vertex model 
is zero [s^. All points in the disordered phase have this line as their RG stable fixed point 
(all points on the line are equivalent). This is the infinite temperature limit of the classical 2D 
eight-vertex model. In this regime the expectation value of the polarization operator {t{A)t{B)) 
is zero (and its correlation vanishes at a length scale of the lattice spacing.). Thus this state 
does not exhibit long-range order, but since the dual variable {ii{A) fj,{B)) is non- vanishing, it 
exhibits topological order. Quantum mechanically, the wave function is the equal-amplitude 
superposition of all configurations of arrows consistent with the eight-vertex restrictions. This 
state is deconfined since in this limit the Hamiltonian reduces to the flip term (plaquette). It 
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is well known 31] that in this regime the Wilson loop operator has a perimeter law behavior 
and that the energy of two static sources is finite and independent of their separation. The low 
energy (i. e. long-distance) sector of this phase describes a topological Z2 deconfined theory, 
equivalent to the Kitaev point. 

3. Critical Behavior. 

The 2D quantum theory has lines of critical points = d? + 2 (and d? = (? + 2). There are 
also lines of critical points for d = and < < 2 (and c = and < < 2). All of these 
critical lines are (up to duality transformations) equivalent to the six-vertex model. Notice 
that, since the energy of this quantum state is exactly zero (by construction), the ground-state 
energy does not have singularities at the phase transitions. This is of course a peculiarity (or 
rather a pathology) of this model and it certainly non-generic: any perturbation leading to a 
non-vanishing energy will lead to singularities in the energy of the ground state. 

In section UlIBl we discussed how the ground-state correlators of the (multi) critical square-lattice 
quantum dimer model were the same as those of the quantum Lifshitz theory, Eq. ()2.4() . at = 2tt 
|l4ll53| . The same equivalence holds for the quantum eight- vertex model for all values of k. Moving 
along the critical line results in changing k as given in (|4.9l4.1()l) . Kadanoff W^I^^ showed that 
the critical behavior of the classical eight-vertex model can be mapped exactly to the critical behavior 
of the two-dimensional Gaussian model, the free boson described by Eq. 1)2.12(1 and Eq. (|2.13() . The 
equal-time correlators in the quantum Lifshitz model summarized in Appendix^are therefore identical 
to those of the quantum eight-vertex model on its critical lines. Moreover, when d = and (? < 2, 
the eight-vertex model becomes the six-vertex model, which has a height description. Thus the same 
heuristic arguments applied in the last section can be applied here, again implying that the quantum 
six-vertex model and the quantum Lifshitz theory are in the same universality class. In particular, the 
dynamics of the quantum eight-vertex model along its critical lines will obey scaling with a dynamic 
exponent z = 2. 

To compete the identification, we need to express the coupling constant k of the quantum Lifshitz 
model in terms of the eight-vertex parameters c and d. This is easiest to do by computing the dimension 
X of the "energy" operator, which when added to the action of (classical) critical theory, moves it away 
from criticality. Combining universality with Baxter's exact result ()4.1|) for the correlation-length 
exponent gives x = 2(1 — /u/vr). We have normalized the boson 93 in Appendix [SI so that the energy 
operator is given by cos(2(/9), which has dimension x = 1/{2ttk). Note that at the free-fermion point 
= 27r, the energy operator has dimension 1, the dimension of a 2D fermion mass term. Combining 
the two and using a simple trigonometric identity yields 



and by duality 



= 8 cot"^ (cd) . for|c^-d^| = 2 (4.9) 



= 8cot"M a/^ - 1 forc^<2,d = 0. (4.10) 




At the Kosterlitz-Thouless transition point = 2, d = where the two critical lines meet, the 
dimension x = as expected, and both formulas give = 47r here. 

This equivalence to a scalar field theory can be extended to the scaling region near the critical 
lines, by repeating the arguments of section IIII CI The perturbing operator is cos(2(/?); symmetry 
forbids cos{(p) from being added to the action of the classical theory. Thus near to the critical lines the 
effective field theory will be of the form (|3.8() and Eq. (|2.9|) : these field theories will have the correct 
critical exponents. 
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C. Operators 

Equal-time correlation functions of operators which are diagonal in the arrow representation of the 
eight-vertex model are identical to correlation functions of the classical eight-vertex model. Unfortu- 
nately, these correlation functions can be computed exactly only along the free-fermion lines, by using 
Pfaffian techniques. Likewise correlation functions of off-diagonal operators {e.g. flip operators) can 
only be determined at some special points deep in the ordered and disordered phases by using directly 
the quantum eight-vertex Hamiltonian. However, at or near the critical lines, the asymptotic behavior 
of the correlators can be found by using universality. The operators of interest in the eight-vertex 
model can all be expressed in terms of the charge and vortex operators of the Gaussian model [95| . 

The operators of interest in the quantum eight-vertex model can be readily identified in the quantum 
Lifshitz model. Thus we have a full identification of a 2 + 1-dimensional quantum critical theory. As 
a direct consequence of being able to identify these operators we can also study how these operators 
perturb the critical theory, at least within a renormalization group argument along the lines of Ref. 
jSJ]; the analog in the quantum dimer model is our analysis of H^^ for w small at the end of section 2. 

To identify the eight- vertex operators in the bosonic language as in Ref. |95l |. we first study the 
six-vertex line, where there is a height description. There the height is local in terms of spin variables, 
so we expect the boson 93 to be local as well. Standard bosonization techniques [o^ show that only 
products of spin operators from both Ising models can be written in terms of the boson. One finds that 
cos{(p) ~ t{A)t{B). Note that this is consistent with the identification of cos(2(/3) = 2cos^(93) — 1 with 
the perturbing operator; the usual Ising operator product gives the fusion rule t{A)t{A) ~ 1 + £{A), 
where £{A) is the energy operator. 

To study the operators in more depth, we utilize a 2+ 1-dimensional Z2 gauge theory of the quantum 
eight- vertex model derived in Appendix IC 11 The degrees of freedom of this gauge theory reside on 
the links of the square lattice and have a one-to-one correspondence to arrow configurations of the 
eight-vertex model. The Hamiltonian of the Z2 gauge theory is 

-f^qSv = Hy + i^flip (4.11) 

where -fTmp has the form of a sum over plaquettes of Z2 flip operators^: 

Hfiip = -^(t3(x) +ei) al{x) cT?(f + 62) (4.12) 

X 

The potential energy terms are combinations of operators which project onto the allowed select a, b, 
c and d vertices, and assign the weights a, b, c and d to different contributions to the wave function. 
The explicit form of the potential energy terms Hy is given in Eq. (|C6|) . As in all Z2 gauge theories 
Isills^ , the physical states of this theory satisfy the constraint of gauge invariance ( "Gauss Law" ) 

G{x) = al{x) al{x) ct|(x - ei) al{x - 62) = 1 Vf (4.13) 

which in this context simply expresses the restrictions on the configurations allowed in the eight-vertex 



model. The operator G{x) is the generator of local time-independent gauge transformations |3lll57 |. 

In Appendix IC 21 we derive the dual theory of the quantum eight-vertex Hamiltonian of Eq. 1)4. 3() . 
The degrees of freedom of the dual theory are defined on the sites of the dual square lattice. The 
Hamiltonian is 



Throughout this section and in Appendix |n| we use the standard notation used in lattice gauge theories E^. The 
(t's are Pauli matrices; the superscript is a Pauli matrix label and the subscript indicates the spatial direction, 1 — 
horizontal and 2 = vertical. Notice that the notation used by Kitaev 42] is somewhat different. Please see Appendix 
im and figures |S] and |U] for details of the notation that we use here. 
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where ^^v.duai is given in Eq. (|C6|) and Eq. HC15|) . We will use the both the gauge theory and its dual 
to investigate these phases. 

Let us discuss briefly the observables of this model, how they behave in the ordered (confined) and 
disordered (deconfined) phases and what critical behavior they exhibit at the phase transition lines. 
As we note in Appendix IC 21 the quantum eight- vertex model retains the two-sublattice structure of 
the classical model. The Z2 x Z2 symmetry of the classical model also survives in the quantum theory. 

1. Order and disorder operators: We denote by t{A) and t{B) the (order) operator for sublattice 
A and B respectively. In the gauge theory this operator is a disorder or kink operator [sj] ; in this 
2D gauge theory this is just the flux (monopole) operator of Refs. |57ll97| . or the vison operator 
of Refs. 0111]. The order operators are just the order parameters of the sublattice Ising models. 
As such they acquire an expectation value in the ordered phase (where the staggered polarization 
has an expectation value as well). In the disordered phase their equal-time correlation functions 
decay exponentially with distance (as do their connected equal-time correlators in the ordered 
phase). The precise behavior of this correlation function is only known on the decoupling curves 
where they reduce to the correlation functions of the 2D classical Ising model, and on the critical 
lines! The scaling dimension of this spin field (which is known as the "twist" field in conformal 
field theory [48il96j ) . is 1/8 not only at the decoupling point (as is well known from the 2D Ising 
model) but along the entire phase boundary. We note however, that the spin fields can not be 
represented directly in terms of the boson, but only in its orbifold 99]. 

Similarly, we will denote by /u(^) and fJ,{B) the "frustration" (or "fractional domain wall") or 
disorder operator 32] for sublattices A and B respectively. In the 2D classical theory, this is 
the dual of the order operators and in conformal field theory it is also represented by a twist 
field. It has the same scaling dimension as the spin field. On the other hand, in the gauge 
theory, the disorder operators correspond to a local violation of the Gauss-Law constraint, Eq. 
()4.13() . by imposing that G{x) = — 1 just at site X, I.e. cL static source (i.e. an "electrode", or 
Polyakov loop) at site x. This operator is analogous to the holon (or "monomer" ) operator in the 
Rokhsar-Kivelson quantum dimer model. This expectation value of products of two operators 
of this type serves as a test for confinement. The disorder operators have a non-zero expectation 
value in the quantum disordered phase. A simple minded calculation shows that, in in the 
ordered phase, the ground state with two such defects (or sources) separated at a distance R, 
has a non-zero energy U{R), which grows linearly with R, U{R) ~ aR, consistent with the fact 
that the defects create a fractional domain wall in an ordered state. We expect that the string 
tension a vanishes as the critical lines are approached, with a behavior dictated by that of the 
gap scale, i.e., Thus, the excitations created by the disordered operator are confined in the 
ordered phase^. In contrast, in the quantum disordered phase, the energy of a pair of defects 
saturates to a finite value if ^ Hence, the disordered phase is deconfined. 

It is also interesting to ask what is the energy cost of a set of vortex configurations on the 
quantum critical lines. This can be done in a number of ways. We note here that using path- 
integrals the change of the ground state energy can be found by computing the path integral 
in a background of vortices, normalized by the path-integral without vortices. Consider the 
simple (an general) case of two vortices of magnetic charges itm separated a distance R. Let 
is denote this amplitude by W{R,T), where T is the (infinite) time-span of the system. In the 
dual gauge theory this is the same as the computation of two static Wilson (or Polyakov) loops, 
corresponding to two static sources with electric charges itm at a distance R from each other. 



^ The holon operator of the quantum dimer model has the same behavior in both the columnar and staggered phases. 
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The energy cost of these defects is 

U{R) = lim (^-^ In W{R,T)^ (4.15) 

We can compute U (i?) by using the path-integral for the quantum Lifshitz theory of Eq. (|2.12|) 
in a background of vortices, which amounts to modifying the action by the minimal coupling 
shift 

(VV)^ ^ (V • {V4> - A)f (4.16) 
By explicit calculation we find that 

U{R) = (4.17) 

In other terms, the interaction energy of vortices is zero! A simple way to understand this 
results is to recall that, in imaginary time, the quantum Lifshitz model is a theory of a smectic 
liquid crystal. The vortex state is just a configuration of screw dislocations running along the 
z axis (imaginary time). Since the smectic is isotropic and has no resistance to shear in the 
xy plane, these defects do not cost any energy. From a quantum mechanical point of view this 
result is consistent with the known fact that at the RK point of the quantum dimer model, the 
interaction energy of a pair of monomers is zero^". In gauge theory language this means that 
at this critical point vortices are completely free.^^ In contras t, at quantum critical points of 
Lorentz invariant gauge theories, which have z = 1, it is known '10(j\ that at their critical point 
the effective interaction has the universal law U{R) ~ 1/R, in all dimensions. 

2. Polarization: The polarization operator, the local arrow configuration on a link, is the natural 
order parameter. In terms of the order operators it reads P = t(A)t(B). As we saw above 
in the classical "antiferroelectric" ordered phase, P is a natural order parameter, although 
appropriately-staggered expectation values of t{A) and t(B) also do not vanish in this ordered 
phase. The behavior of the polarization operator on the critical lines of the classical model was 
studied by Kadanoff and Brown [o^, who found that its scaling dimension is Ap = 1/8itk, with 
K given by Eq. (|4.9j) . From the results of Appendix^ we see that we can identify the polarization 
operator with the boson operator P ~ cos((^) in the quantum Lifshitz theory, which has the 
same scaling dimension. 

3. Mass term: In the classical theory the mass term, or "energy density", is the product of two 
(dual) spin variables r on nearest neighboring sites on the same sublattice. We will denote them 
by £{A) and £{B) respectively. Although in the quantum theory these operators no longer 
correspond to an energy density, we will define £ = {£{A) + £{B)) /\/2, as usual, and refer to 
this operator as to the mass term since it drives the theory away from criticality. In Ref. 

it is found that the scaling dimension of this operator is Af = l/27rK = 4Ap, which led to the 
identification £ ~ cos(2(^). If we regard this operator as acting on the critical line along the 
d = axis, i.e. the 6- vertex model, we see that this is a relevant operator for the quantum 
Lifshitz model as well. Grinstein's RG arguments [H^] show that this perturbation drives the 
system into a massive phase with a residual unbroken Z2 symmetry. This is just the disordered 
phase of the quantum eight-vertex model, and this perturbation drives the system towards the 
Kitaev point. A similar analysis works on the = -|- 2 critical line. 



We thank Shivaji Sondhi for this remark. 

This result, Eq. 14.171 holds strictly only at the quantum Lifshitz fixed point (or rather line) fixed point, Eq. (12.121 . 
However, marginally irrelevant operators such as (V(^)*, discussed extensively in Refs. ,,75..,76l . give rise to corrections 
to scaling which change the behavior of U{R) from Eg. 14. 171 to a (ln_R)/_R^ law||75|. We thank T. Senthil for pointing 
this out to us. 
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4. Wilson loops: As it is well known gauge theories, the natural test for confinement is the behavior 
of the Wilson loop operator = Ylg^^ 'Zfl^)' "^^ere {£} is a set of links that belong to a closed 
path 7 on the (direct) lattice (see Ref. [53] and references therein). Since these operators are 
off-diagonal in the arrow representation, i.e. in the representation in which is diagonal, they 
do not have an analog in the classical 2D eight-vertex model. A simple argument [sj] shows 
that in the ordered state, which for c — > oo is essentially an eigenstate of a^, the expectation 
value of the Wilson loop is zero and that first non-zero contribution comes from acting A^(7) 
times with the flip operator, where A^(7) is the number of plaquettes enclosed inside the path 7, 
and that the behavior thus found has the form exp(— constant N('j)), i.e. an area law. Thus the 
ordered state is a confining phase. Conversely, deep in the disordered phase, the ground state 
is essentially an eigenstate of the flip operator and we get a perimeter law, i.e. deconfinement. 

5. Spinors: The classical 2D eight-vertex model has two Majorana fermion operators ^'{A) and 
tp{B), constructed as usual as products of order and disorder operators [s^l, ip ~ TfJ-- We can 
also define an analog of the spinors in the quantum theory (see the discussion in Appendix ^ 
on charge and vortex operators in the 2 + 1-dimensional quantum Lifshitz model). The scaling 
dimension of the spinors is [o^ = (27rK + l/27rK) /4, a result familiar from the Luttinger 
model. 

6. Other operators. The classical 2D eight-vertex model contains a marginal operator E = 
£{A)£{B), with scaling dimension A^; = 2, which is responsible for the lines of fixed points 
with varying critical exponents. The quantum eight-vertex model discussed here has a similar 
behavior. In the quantum Lifshitz model the marginal operator is^^ E ~ (V^V') . Another 
operator of interest is the crossover operator {S{A) —£{B)) /\/2 which breaks the sublattice 
symmetry and has scaling dimension 2ttk, The operator product expansion of two crossover 
operators generates an operator symmetric under sublattice exchange, and is identified with a 
charge-two vortex operator. The scaling dimension of this operator is Sttk and it is becomes 
marginal at k = l/47r, driving the Kosterlitz-Thouless transition along the d = line.^^ Other 
interesting operators are the "mixed" operators P = t{A)h{B) and P* = fi{A)T{B), both with 
scaling dimension Ap = 7rK/2. Hence P and P* are vortex operators as well; these are non-local 
with respect to P. They have been identified as the holon creation operators in the context of 
the quantum dimer model 



Let us end this section by noting that there is a simpler quantum eight-vertex model that one can 
write down whose behavior is quite different from the one we have discussed here. Basically, this model 
will give weights to all the vertices, according to their type. In addition, there will be a flip term. In 
its dual form the model is 



where the projection operators are defined in Eq. HC15|) . The conventional Baxter weights of the 



classical model {i.e. when the flip term is absent) are a = e *s^, etc. 



As we noted above, on symmetry grounds a lattice model may generate additional opera tors w hich break rotational 
invariance, absent at this RK point, which may drive the transitions to be first order, see l75l76|l . 

Notice that the charge-one vortex operator is not allowed in this case, which is the operator driving the KT transition 
in the 2D classical XY model |101iil0ail03Ml04nl . 



D. Another quantum eight-vertex model 




(4.18) 
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We first note that this model is not of the Rokhsar-Kivelson type (except when Sa = Sf, = Sc = e^, 
which is a dual version of the Kitaev model). The ground state wave function does not look like the 
statistical weights of the classical two-dimensional model, i.e. it is not a linear superposition of states 
weighted according to the vertices present. Thus, strictly speaking we can not use the machinery of 
this paper to study the ground state of this model. We do expect, however, that this model will have 
the same ordered phases, just like the ferro-electric and anti-ferro-electric phases of the eight-vertex 
model, as well as the same confinement-deconfinement properties discussed above. In addition, we 
expect that these ordered phases have the same spectrum and general properties as the ones we found 
in the Rokhsar-Kivelson type model we studied above. Similarly, in the regime where the flip term 
dominates, the ground state of this model is a uniform quantum disordered state, with the same 
properties as those of the topological phase of the model studied here. 

However, the quantum critical properties of Eq. (|4.18|) are completely different from the quantum 
eight-vertex model discussed in the rest of this section. They instead have a number of features in 
common with well-known three-dimensional classical models. By setting Sa + £b = + ^d, this model 
reduces to two interpenetrating two-dimensional Ising models in a transverse field, as the four-body 
interactions cancel. Thus along these curves this model has second-order phase transitions in the 
universality class of the classical three-dimensional Ising model, as opposed to the quantum eight- 
vertex model discussed above. The latter has z = 2, while the model of Eq. (|4.18() therefore has 
dynamical exponent z = 1 and is Lorentz invariant in the critical regime. It is easy to see that this 
behavior extends beyond the decoupling point. Indeed, the universality class of the 3D classical Ising 
model is controlled by the Wilson-Fisher fixed point and is accessible by the 4— e expansion. Away from 
the decoupling point we have then two 3D Ising models which are coupled through their energy density. 
This theory is in the universality class of a two-component real field (ipi,ip2), whose U{1) symmetry 
is explicitly broken down to a discrete Z2 x Z2 symmetry by quartic operators of the form (Pi(p2- This 
is a cubic symmetry-breaking perturbation and it is well known to be perturbatively relevant at the 
Wilson-Fisher fixed point of the two decoupled Ising models, but irrelevant at the Wilson-Fisher fixed 
point for the classical XY model^^. It is well known from perturbative 4 — e RG studies that the 
transition near the decoupling point that the transition may also become a fluctuation - indu ced first 
order transition depending on the sign of the effective coupling (see for instance Refs. |3 ll05| ). 

V. NON-ABELIAN TOPOLOGICAL STATES 

In this section we broaden our scope and discuss field theories with continuous non-abelian sym- 
metries. One can of course obtain a quantum critical point with non-abelian symmetries by taking 
copies of the quantum Lifshitz Hamiltonian. By choosing the charge lattice of the vertex operators 
carefully one can describe a theory with a non-abelian symmetry, say SU{N). However the equivalent 
conformal field theories thus obtained always have integer central charge and are at level 1. Although 
critical, such theories only support excitations with abelian statistics. One may try to attempt to 
generalize this result to models by writing down a theory with more structure, such as one related to a 
Wess-Zumino-Witten model. However, we will find that, contrary to naive expectation, such theories 
will turn out not to be critical but, instead, in a phase with a finite gap, a topological phase, just 
like the fractional quantum Hall effect. We will see that it is not at all clear how to find a conformal 
quantum critical point with a non-trivial non-abelian symmetry. 

Nevertheless it is possible to use the approach of the previous sections to discuss topological non- 
abelian phases, even if any associated non-trivial quantum critical points are not yet known or do not 
exist. By analogy with the theory of fractional quantum Hall states, as well as from general results 



This is a classic example of a dangerously-irrelevant operator, since it removes the would-be Goldstone boson from the 
broken-symmetry phase. 
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in Chern-Simons gauge theory, one expects that a time-reversal-breaking topological phase in 2 -|- 1- 
dimensions should have an effective field theory description in terms of a Chern-Simons gauge theory 
at some level k. On closed manifolds this is a topological field theory in the sense that the partition 
function is independent of the metric of 2 -|- 1-dimensional space-time, and that the expectation values 
of its gauge-invariant observables, the Wilson loops, depend only on the topology of the loops such as 
the knot invariants |2^. 

By analogy with the approach pursued in the previous sections, we will discuss the properties of 
the wave functionals of this theory. This will lead us to consider first the strong-coupling limit of 
Yang-Mills theory with a Chern-Simons term, and later the topological sector of this theory. We will 
find the ground-state wave functional of this theory, and show how it is related to a two-dimensional 
Wess-Zumino-Witten (WZW) model for a specific choice of polarization [i.e., the choice of canonically 
conjugate variables) which is natural in the strong coupling limit of the Yang-Mills Chern-Simons the- 
ory. Even though the WZW model is a conformal field theory with algebraically-decaying correlators, 
the gauge-invariant ground-state correlators in the two-dimensional quantum theory are exponentially 
decaying. Topological invariance will only be attained by correlators of gauge-invariant operators, 
Wilson and Polyakov loops, and only at long distances. The excitations obey non-abelian statistics: 
when bringing particles around one another, not only does one pick up phases, but the order in which 
they are exchanged affects the final state. 

We will also study a doubled Chern-Simons theory which, contrary to the fractional quantum Hall 
effect, is time-reversal invariant. Using by now standard methods, we will construct the ground-state 
wave functionals of this doubled theory for a gauge group G at level fc, and use its connection with a 
gauged WZW model to understand its properties. We note that the doubled theory is often used to 
circumvent technical problems involved in the use of the formally anomalous wave functional of the 
undoubled theory: see the discussion in the subsections IV Bl and IV CI The topological properties of 
the doubled theory have been recently discussed in depth by Freedman et al. [50|. 




Let us attempt first to construct a non-abelian version of the quantum Lifshitz theory, which one 
might hope would still be a quantum critical theory. We will follow the same procedure as in the abelian 
case. Since there are many two-dimensional Euclidean critical points with non-abelian symmetry, the 
first thing to try is the above procedure for making a two-dimensional quantum theory from a two- 
dimensional classical theory. This procedure is indeed fairly general. One can easily construct a set of 
a appropriate projectors Qi in a lattice model without local constraints. However, the caveat "without 
constraints" is quite important: many interesting classical lattice models have local constraints like 
in the dimer model and in the six- and eight-vertex models. This means that one must find a flip 
operator which respects the constraints. There is no guarantee that there is a set of flips which are 
both local and ergodic. This is particularly apparent in models where the two-dimensional degrees 
of freedom are closed loops. At a critical point, arbitrarily-long loops are an important part of the 
configuration space, so one must construct some sort of flip which still acts non-trivially but locally 
on these loops. In the continuum limit, it is easy to imagine such flips, but it is not always obvious 
how to make them act consistently on the lattice. Nevertheless, the Rokhsar-Kivelson Hamiltonian 
and the quantum eight-vertex model provide examples of models where one can solve this problem. 

For fleld theories, the issue is similar. If a theory has an action S2d and one ignores constraints such 
as gauge invariance, one can easily flnd projection operators Q{x) which annihilate the state weighted 
by e~'^2d^ Por a single fleld we have 




A. Field theories with continuous non-abelian symmetries 




^S2d 




and the Hamiltonian is deflned via Eq. (|2.9|) . This Q is hardly unique: one can multiply it on the left 
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by any operator without changing the ground state. It is important to note that this construction is 
not always easy to implement. Many interesting two-dimensional field theories do not have a simple-to- 
describe action: they are usually defined instead in terms of representations of an infinite-dimensional 
symmetry algebra such as the Virasoro algebra, arising from conformal invariance [i^. For example, 
the action of a G/H coset model [e.g. a conformal minimal model) is written as a gauged WZW model 
such as we will discuss below; one must then worry about the constraints of gauge invariance. 

A point we wish to emphasize, however, is that even if one succeeds in constructing such a Hamilto- 
nian, the 2-|-l-dimensional theory may be qualitatively very different from the classical two-dimensional 
theory. For example, the classical theory may be critical and have algebraically-decaying correlators, 
but the quantum theory may have exponentially-decaying correlators in the ground state. The reason 
is that quantum-ground-state correlators are weighted by l^'oP, not ^o- For the models discussed 
above, this has effects which are easy to take account of: it is the reason that the axes in figure [H] are 
labeled by (? and instead of c and d. If however ^'o has a complex part, the effects can be much 
more dramatic. 

Let us first discuss the simplest two-dimensional classical field theories with a non-trivial non- 
abelian symmetry. The G principal chiral model is written in terms of a unitary matrix field g taking 
values in some simple Lie group G, with action 

SPCM = j d'xTr [g'^d'^gg'^dag] . (5.2) 

This action has a global G/, x Gr symmetry under g LgR) , where L and R are elements of separate 
groups dubbed Gl and Gr. This model looks critical: the coupling constant u is naively dimensionless. 
However, in two dimensions, the one-loop beta function is proportional to the curvature of the target- 
space manifold G. Simple Lie groups G are curved: for example, SU{2) as a manifold is isomorphic 
to a three-sphere. Moreover, one finds that the trivial fixed point at u = is unstable, and the model 
has a finite correlation length proportional to an exponential in \/u. Correlators decay exponentially, 
not algebraically. The phenomenon where a naively dimensionless coupling becomes dimensionful 
due to loop corrections is known in the particle physics literature as dimensional transmutation; it 
is familiar in condensed-matter physics in the Kondo problem and the Hubbard model at half-filling. 
Another reason a gap should appear is that the Mermin-Wagner-Coleman theorem does not allow 
Goldstone bosons in two Euclidean dimensions. One would obtain SpcM in a theory as the low-energy 
limit of a fermionic theory where one attempts to spontaneously break the chiral symmetry Gl x Gr 
to its diagonal Gd subgroup. The would-be Goldstone bosons would take values on the manifold 
Gl X Gpi/Go — G. This satisfies the theorem by having the low-energy excitations get a mass and 
restoring the symmetry to the full Gl x Gr. 

A critical theory with non-abelian symmetry G does occur when an extra term, the Wess-Zumino 
term, is added to Srcm- This term is easiest to write in a three-dimensional space M which has 
the two-dimensional space of interest as a boundary; we assume the two-dimensional space has no 
boundary. It is 

r(ff) = / d'^ e^^'^Tr [g-'d^g g-'d^g g-'d,g] (5.3) 

One finds that different ways of extending space to three dimensions result only in changing T{g) by 
2tt times an intege r. T hus we can add ikT{g) to Srcm for any integer k. This model has a critical 
point at u = 'i/\k\ |l06| |. This theory, with action 

I{g) = SpcAiig, u = l/\k\) - ikT{g) , (5.4) 

has conformal invariance and is known as the G^ Wess-Zumino- Witten (WZW) model. A great deal 
of information is known about the WZW model; the most important for our current purposes is that 
the correlators decay algebraically, as they must at a critical point. 
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Since the WZW model has an explicit action, it is easy to use the trick in Eq. (|5.1|) to find a 
two-dimensional Hamiltonian which has a ground-state wave functional with amplitudes given by the 
WZW action, i.e., 

^o[9]=e-'^'l (5.5) 

The WZW model is critical, so one might guess that the quantum version with this wave function 
would also be critical, in the fashion of the theories in the previous sections. However, it is not. The 
reason is that the Wess-Zumino term has an i in front of it, so 

|^0[5]|2 = ^-2SpcM{g,u=l/\k\) ^^Q^ 

Thus equal-time correlators in the ground state are weighted with the action of a theory with a finite 
correlation length, the principal chiral model. These correlators decay exponentially. 

We do note that when the level k = 1, one can use an alternate method to find the 2d WZW 
correlators in the ground-state of a quantum theory. All the fields of SU {N)i can be written in terms 
of N — 1 free bosons |4^ . The combination of the corresponding quantum Lifshitz theories will then 
have the SU{N)i correlators in its ground state. Likewise, 0{N)i and SU{2)2 can be written in terms 
of free Majorana fermions. Since the different bosons or fermions here do not interact, the physics is 
the same as that discussed earlier; in particular the model is essentially free, and the excitations are 
abelian. 



B. Yang-Mills, Chern-Simons, and Wess-Zumino- Witten 

Since it is not possible to have a wave functional with [^'[(jjp = e~^^^\ we will need to work harder 
to find how the physics of the WZW model can arise in 2 + 1 dimensions. A connection between 
the WZW mod el an d the wave functions of 2 + 1-dimensional Chern-Simons gauge theory has long 
been known Chern-Simons theory without any matter fields is a topological field theory: 

the physical states are Wilson and Polyakov loops, and their correlation functions do not depend on 
distance but only on topological properties of the loops. These correlators can be expressed in terms 
of the Verlinde numbers, which also describe the (chiral) fusion rules of the WZW model ^lo3 |. These 
topological field theories have been discussed in detail in many places. One discussion of issues closely 
related to those of interest here is the recent work of Freedman et al. (s^l ; we will discuss the relation 
of our results to this paper in the next subsection. 

In this subsection we will describe a wave functional which after squaring and averaging over fields 
indeed yields the WZW partition function. We will show that this wave functional is the ground state 
of a 2+1-dimensional gauge theory. This all sounds like exactly what we want for a quantum critical 
point, but in a gauge theory, physical observables must be gauge invariant. Gauge-invariant states 
will turn out to have a gap, and gauge-invariant correlators in the ground state will be exponentially 
decaying. Introducing a gauge field therefore will not yield a quantum critical point with non-abelian 
symmetry. It will yield, however, something just as interesting: an explicit Hamiltonian for a theory 
in a topological phase. 

Before discussing construction of these wave functionals, which was done in great detail by Witten 
[s^ l and whose construction we will follow here, it is worthwhile to discuss first its physical meaning 
and to emphasize some well-known facts of Chern-Simons gauge theories 26|. The ground state wave 
function of any field theory can be viewed as the quantum mechanical amplitude of some arbitrary 
state into the vacuum (or ground state), i.e. what we normally call the vacuum state in a given 
representation. As such this amplitude is the functional integral of the quantum field theory on an 
open manifold bounded by the initial time surface. However, the partition function of a topological 
gauge theory such as Chern-Simons is only gauge invariant on a manifold without boundary. Thus, 
the ground state wave function has a gauge anomaly. Nevertheless, if the wave function is used to 



30 



compute expectation values of gauge-invariant observables, the result is gauge invariant. The reason 
is that the computation of expectation values (as well as all inner products) involves the conjugate 
wave function, i.e. the amplitude to evolve from the vacuum state in the remote past into the chosen 
state at the fixed time surface. Consequently, the gauge anomaly cancels in the computation of the 
expectation values of gauge-invariant operators. In this way, at a formal level, the computation of 
expectation values leads to a formally "doubled" theory even though the number of degrees of freedom 
has not changed. In subsection IV CI we will discuss a physically doubled theory which has a formal 
relation with what we do in this Subsection. There is an extensive literature on the technical issues 
involved in this problem: the relevant discussion for the analysis done in this paper can be found in 
Refs. |52..109...110.111j. 

We should note here that the gauge anomaly of the wave function is formally (mathematically) 
analogous to the gauge anomaly of Chern-Simons theory on manifold with a spatial edge. In the latter, 
the anomaly is physical: to cancel the anomaly the gauge-invariant theory must include physical degrees 
of freedom residing at this 1 + 1-dimensional boundary. In the incompressible fractional quantum Hall 
fluid, this results in physical edge states 1121 113| . 



Another important fact is that, since Chern-Simons gauge theory is a topological field theory, the 
partition function on a closed space-time manifold is independent of the metric. However, the boundary 
of the manifold and the choice of polarization break the general coordinate invariance of Chern-Simons 
theory. In particul ar, the choice of holomorphic polarizatio n in duces a conformal structure in the wave 
functional 2^ 107l | which is absent in other polarizations 114 |. We will see that the conformal wave 



function arises naturally in a specific theory, the strong-coupling limit of the Yang-Mills Chern-Simons 
theory. 

The wave functional of interest was discussed in depth by Witten It involves the WZW field 
g coupled to a gauge field Ai taking values in the Lie algebra of G. It is most convenient to write the 
gauge field in complex coordinates Az = {Ai — iA2)/2 and A-^ = {Ai + iA2)/2. When the symmetry 
Gr of the WZW model is gauged, the closest thing to a gauge-invariant action is 

I{g,A) = I{g) + d^zTv [2A-g-^d,g - A^A,] (5.7) 



This action is not gauge invariant; the subgroup Gr of Gl x Gr is anomalous [5^. Indeed, the fields 
under gauge transformations U{z,z) in Gr as 

g^gU, Ai^U-^AU + U~^diU (5.8) 

so that the action I{g, A) transforms as 

lig, A) ^ I{g, A) + ^J Tr [A.UchU-^ - A^Ud.U-^] - ikT{U). (5.9) 

Note, however, that the variation of I{g, A) under gauge transformations is independent of g. If we 
were to simply square a wave functional proportional to e~^^^'^\ the WZW term in the action would 
cancel and we would be back a gauged principal chiral model. Instead, we define a wave functional 
depending only on A by doing the path integral over g: 

^[A] = j [Dg\e-^^9A)_ (5.10) 

This wave functional does yield the full WZW partition function after squaring and integrating 
over A [s^]. We define this integrated as 

1^1'^^ l[DA]MA]\\ (5.11) 
vol G J 
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where we have divided the measure by the volume of the gauge group G because it follows from Eq. 
(|5.9() that even though ^[A\ is not gauge- invariant, ^'[A]^'[j4] is. Substituting the definition Eq. ()5.1U() 
gives 

= ^ l[DA][Dg][Dh-']e-^^9,A)-{I{h-\A)Y 
vol G J 

I{g,A) + {I{h-\A)r = I{g)+I{h) + ^ j ^f^Ti [A-g-^d,g - A.d^h ■ h-^ - A- A,] 

Note that we have defined A so that (A^)^ = —A^, i.e. covariant derivatives have no i in them. We 
have also used the easily-proven fact that I{h) = I*{h^^). Since the integral over A is Gaussian, it 
can easily be done, giving 

\^\^ = ^^J[Dg][Dh-']exp(^-I{g)-Iih) + ^ j <fzTv[g-^d,gchh ■ h'^]^ . 

This can be simplified by using the Polyakov-Wiegmann formula |ll5| | 

I{gh) = I{g) + I{h) j Tr [g'^d^g d^h ■ h'^] . (5.12) 

Thus the integrand depends only on the product / = gh, so we can change variables to f , which 
cancels the volume of the gauge group. This yields the final expression for the normalization 52!| 

|^|2 = [ [L>/]e-^(^) . (5.13) 



Thus, remarkably, the integrated square of the wave functional ^[A] ends up giving the full WZW 
path integral, including the imaginary piece: because of the path integral, the right-hand-side of Eq. 
()5.13|) is positive and real as it must be. 

To reemphasize a point made earlier, 1)5. 13(1 does not mean we have now found a 2 + 1-dimensional 
quantum critical point whose equal-time correlators are those of the WZW model. We must first find a 
Hamiltonian which has ^[A] in Eq. ()5.1U|) as its ground state. Once having found that, the physically- 
relevant operators are only those satisfying the proper gauge-invariance properties and regularization. 
After having done so, we will end up seeing that despite sharing some key properties with the WZW 
model, the 2 + 1-dimensional model is gapped and has exponentially-decaying correlators. 

To find such a Hamiltonian, let us go back to the definitions Eq. 1)5. 7|) and Eq. (|5.1U|) of the wave 
functional Because the only dependence on A^ is through the quadratic term, we have 

A^-^ ^A] = 0. (5.14) 



6Az 47r 

With a little more work l52ll, one also can show that 

where we have define the covariant derivatives via Dig = dig — gAi and DiAj = diAj + [j4j,^j]. The 
field strength is defined as F-^z = dzA^ — dzAz + [viz, Az\. 

We would therefore like to find a Hamiltonian with a ground-state wave function satisfying Eq. 
(|5.14|) and Eq. 1)5. 15|) . As noted in Ref. [s^ ]. these equations arise in the canonical quantization of 
Chern-Simons theory |l07l ]. so this suggests we look there. The precise Hamilton ian turns out to be 
given by the strong-coupling limit of Yang-Mills theory with a Chern-Simons term |ll(]| | . In the action. 
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we have a gauge field taking values in the Lie algebra of G; here /i = 0, 1,2 and A^ depends on 
space and time. The action for the strong-coupling limit of Yang-Mills theory on a three-manifold M 
includes only the electric-field term, namely 



1 



=2^ L.^' 



M 



(5.16) 



This term is not Lorentz-invariant, but does preserve two-dimensional rotational symmetry. The 
Chern-Simons term is 



Tr 



M 



A^di/Aa -\- ^A^AyAa 



(5.17) 



Under gauge transformations U{x) belonging to G, the integrand in Ssc is invariant, but the integrand 
in Scs is not: 



Scs -^Scs + ^ y^ e^-'" Tt [d^ {A^ + UdM-')] + kr{U). 



(5.18) 



For the manifold M a three-sphere, the latter term turn s out to be the winding number of the gauge 
transformation U{x), and is an integer times 27rfc |ll6l |. The Chern-Simons term is gauge-invariant 
if M has no boundary, and k is an integer. If M has a spatial boundary, one must include massless 
chiral fermions on the edge to restore gauge invariance, giving for example the famous edge modes in 
the fractional quantum Hall effect |^. 

The Hamiltonian with 1)5. 1U() as a ground-state wave functional comes from canonically quantizing 
the theory with action 



S = Scs + Ssc, 



(5.19) 



following, for instance, Refs. lOTllllOl lid ]. Since this is a fairly standard computation, we will be brief 
here. The gauge invariance allows us to fix temporal gauge Aq = 0, so that the degrees of freedom are 
the gauge fields Ai and ^2- Their canonical momenta are 



k 

Hj = -irFoi + —eijAj . 



(5.20) 



The Chern-Simons term contributes nothing to the Hamiltonian, because all the terms are first-order 
in time derivatives. The classical Hamiltonian in the strongly-coupled limit is therefore 



H=^ j (fxTi [{Foif] = y t^^a;Tr 



Hi 



k_ 



(5.21) 



Expanding A and F in terms of generators of the Lie algebra of G, we impose the canonical 
commutation relations 



[A'^{x),U'^M]=iSijS'''d'^'\x-y) . 



(5.22) 



One of the effects of the Chern-Simons term is that Ai and A2 do not commute. In the Schrodinger 
picture, Ilj is given by 



^SA-i 



(5.23) 



operating on the wave functionals. 
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The Hamiltonian of this theory in the Schrodinger picture has a very simple form. Because FQi 
includes Ilj, one has the usual ordering ambiguity in the quantum Hamiltonian. We define 



E = ^(F„,+.F„,) = -A.-A^, (5.24) 



SO that 



[E-{x), {E\y)y] = ^5'^Hi^){x-y). (5.25) 

The operators E and E'^ are like annihilation and creation operators. We then normal-order the 
Hamiltonian as in section^ subtracting the vacuum energy to give 

H = e^ [ (fxTr\E^x)E{x)\ . (5.26) 



Thus the Hamiltonian of strongly-coupled Yang-Mills theory with a Chern-Simons term is precisely 
of the form Eq. (|2.9jl . like all the others discussed in this paper. Moreover, the annihilation relation 
£'^'[^] = here is identical to the relation Eq. 1)5.14(1 . We thus have found an explicit Hamiltonian 
whose ground state obeys the first of the two relations satisfied by ^[A] above. 

Finding the Hamiltonian alone does not complete the canonical quantization of gauge theories with 
a Chern-Simons term. Fixing the gauge = still allows time-independent gauge transformations. 
Moreover, there is no time derivative in Eq. (|5.19|) . so Aq should be viewed as a Lagrange multiplier 
which results in a constraint. When we are canonically quantizing the theory in = gauge, this 
constraint is implemented on the wave functions. Specifically, one has 

D-^ + D,^ + -^d^A- - ^^^A^ ^A] = 0. (5.27) 
OA^ OAz 47r 47r / 

This operator is the generator of time-independent gauge transformations on A^ and A-^, so this con- 
dition amounts to requiring that the wave function be invariant under such transformations. Another 
way of viewing this constraint is as requiring that the wave functionals obey the non-abelian version of 
Gauss' Law. Requiring Gauss' Law along with Eq. I|5.14() yields Eq. ()5.15|) . the other desired relation: 
adding Eq. and Eq. ((^1^ together yields Eq. ((^^TTjl . Moreover, the fact (|^ that "^[A] is not 

gauge invariant is precisely the effec t of the fact (|5.18|) that the Chern-Simons action is not gauge 
invariant when there are boundaries jl07j |: one can think of the constant-time slice required to define 
a wave functional as a boundary in space time. 

We have seen that the wave functional Eq. (|5.1U() indeed describes a zero-energy ground state of 
the 2 + 1-dimensional theory with action Eq. ()5.19|) . We now need to understand the corr elators in 
the grou nd state and the excited states. Luckily, the former issue has been studied in Refs. flOflHTTTl . 



nilTTi], and the latter in Refs. [TinlliTil| . We will show in the next subsection that the ground-state 



are those of a topological field theory, and that one can prove there is a gap in the spectrum. 



C. The doubled theory 

It is both convenient and physically important to study the wave function of the "doubled" theory, 
where we have two gauge fields A and B. The three-dimensional action is the sum of an action of the 
form of Eq. (|5.19() for both fields, where the Chern-Simons term for the field A has coefficient k, while 
that for B has coefficient —k. As opposed to Chern-Simons theory with a single field, the doubled 
theory is also invariant under time reversal and parity, if the fields A and B are exchanged under these 
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transformations^^. The two fields are not coupled in the action, so formally the wave functional of the 
doubled theory factorizes: 



X[A,B]^^AMB] (5.28) 

There is also a compelling technical reason to study the doubled theory. We will show that the doubled 
theory amounts to gauging a non-anomalous symmetry. In particular, after integrating out one of the 
gauge fields (a simple Gaussian integration), the wave functional is gauge invariant in the remaining 
field. The theory can be quantized consistently; effectively the regularization (the measure of the path 
integral) couples the two copies. Thus in spite of this factorization of the wave function, the operators 
that create the physical states are made of operators acting on each sector, carefully glued together to 
satisfy the requirements of gauge invariance. 

Canonically quantizing the doubled theory is straightforward, since the action splits into decoupled 
pieces. The ground-state wave functional of the doubled theory can be written as [s^ ] 



X[A,B]^ j [Dg]e~'^3A,B) ^5 29) 

where 

I{g, A, B) = I{g) + cPzTv [2A^g-^d,g + 2B,gd^g-^ - A^A, - B^B, + 2B,gA^g-^] (5.30) 

One can prove that oc ^[j4]^'[i?] indirectly by showing that x[^)^] satisfies both (|5.14|) and 

()5.15|) for A, and the conjugate equations for B. Directly, we prove this by first noting that the 
Polyakov-Wiegmann identity 1)5.12(1 yields 

I{g,A) + {Iih,B)r = I{h-'g,A,B) - ^ J d'^zTr [{gD^g-')ihD,h-^)] . 

where the covariant derivatives are D-g~^ = d^g'^ + A-g~^ and D^h^^ = dzh^^ + B^h^^. In the path 
integral we can split apart the last term by introducing an auxiliary gauge field C: 

^[Ai^B] = J [Dg][Dh-^]e-^'^^'^^-'^^^'''^'>^* 

= I [Dg] [Dh-'] [DC] exp ( - I{h-'g, A, B) 

-A j SzTt [aCzigDzg"^) + aCzihDzh^^) + a^C.C^] ) , 

where a is a (small) coupling constant. The terms linear in C can be removed by redefining g and h. 
Namely, under g ^ g + Sg for small Sg, the action 

I{g,A) ^ I{g,A) + ± ^d^zTV [gD^g-'dzidg g-')] . (5.31) 

Thus if we redefine the fields g and h so that 5g and 5h^^ obey aCz = dz{dgg^^) and aC- = 
—dz{h6h~^), the mixed terms cancel and the field C decouples. The integrand then only depends on 
the combination h~^g, so changing variables to / = h~^g yields B], up to a factor of the volume 
of the gauge group: 

^A]^B]= J [Dg][Dh'^]e-^^''''^"^'^'> =Yol G x[A,B]. (5.32) 



For this reason, this theory has arisen for example in theories of quantum computation an d in supercon ductivity |5Cll5H| . 
and in effective "coset" field theories of the non-abelian fractional quantum Hall states |l2nlll2llll22^ . 
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To understand this wave functional we compute its integrated norm. One can easily do 

one of the two functional integrals over the gauge fields in the same manner as before, yielding [54] 

I [DB] \x[A,Bf = I [Dg]e-'o/c(9,A) (5.33) 
where Ig/g is the action of the G/G gauged WZW model, namely 

Ig/g{9, A) = I{g) + ^Jd^z Tr [A^g'^d^g + A.gd^g-^ - A^A, + A^gA^g^^] . (5.34) 

As opposed to I{g,A) in Eq. (|5.7|1 . this action is gauge invariant, because one is gauging the full 
Gl X Gr symmetry of the WZW model. This means the path integral of the doubled theory is well 
defined and free of anomalies, as well as parity- invariant. 

Let us compare the doubled theory to the undoubled one. For the undoubled theory, we found that 
the effective partition function l^'p, the ground-state wave functional squared and integrated over the 
(one) gauge field, was that of the WZW model. However, constructing excited states and operators 
directly in the undoubled theory appears to be problematic due to the gauge anomaly discussed 
above. Although it would seem that a priori one could not require that operators be gauge invariant 
if the ground-state wave functional itself is not, as we emphasized in Subsection IV Bl from general 
considerations of Chern-Simons theory we know that the physical observables are gauge invariant, 
that only gauge invariant observables must be considered, and that their expectation values are free 
of any anomalies. However, while this is apparent in the path-integral construction of the quantum 
theory, it is not so apparent if one is to use the wave function, i.e. in terms of a chiral Euclidean 
WZW model. As we discussed above, at this level one is led to introduce a formal "doubled" theory 
even for the undoubled theory. Thus, although the wave function itself factorizes the physical states 
cannot be constructed in terms of arbitrary factors from each sector of the doubled theory. Hence, the 
requirement of gauge invariance can spoil the apparent factorization suggested by the wave function 
X[A,B] of Eq. 

The formally correct way to define the correlators in the undoubled theory is in the doubled theory: 
one can always introduce another field B and then integrate it out. The reason for doing this is that 
the doubled theory can be properly regulated, because ()5.33|) is gauge invariant. All the formal 
manipulations done above are well founded, because the path integral can be defined properly. In 
contrast in a truly doubled theory, the gauge-invariant physical observables couple to both fields 
which must then be regarded as genuine degrees of freedom. In contrast, in the undoubled theory, 
the additional degrees of freedom are a formal device used to regulate the theory^^. This means that 
correlators in both theories are those of the G /G gauged W ZW model, which is anomaly-free. These 
correlators have been studied in great detail |l09llll7l|ll8| |. In particular, careful discussions of the 
proper regularization of this theory can be found in these papers. As common in non-abelian gauge 
theories, to properly do these path integrals, one needs to introduce fermionic ghosts. One then finds 
a BRST charge Qbrst obeying {QbrstY — 0? physical states are annihilated by it. We can thus 
consistently demand that states and physical operators be gauge invariant. 

A re lation between the G/G gauged WZW model and topological field theory was conjectured 
in jl09l | and proven in [s^ ]. In particular, it was shown that once we integrate over A as well, the 
resulting partition function |xp is independent of the metric of two-dimensional space. The correlators 
in the ground state therefore are independent of distance, and given by a topological field theory. A 
topological field theory is obtained by studying the states which are annihilated by Qbrst but not 
given by Qbrst acting on something else. In mathematical language, the physical stat es of the 
topological theory are given by the cohomology of Qbrst- It was derived directly in |l23| | that the 



This approach is reminiscent of stone soup. 
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correlators of the G/G topological field theory are given in terms of the Verlinde numbers |l08| |. which 
give the dimensions of conformal blocks in the ordinary WZW model. 

Going to the doubled theory has therefore allowed us to not only avoid technical problems, but also 
to prove that the ground state correlators are those of a topological field theory. We should stress once 
again that the doubled theory used here has twice as many degrees of freedom as the undoubted theory. 
It is a physically distinct time-reversal invariant theory, unlike the undoubled theory which breaks time- 
reversal symmetry. This physical difference is apparent from the construction of its observables: it 
has twice as many "anyons" , which come in time-reversed pairs. On a torus, the Wilson loops of the 
corresponding topological field theory wind around both cycles [HOl. The weights of the loops can be 
defined locally, and one can see how the non-abelian statistics arise |50l |. The results of this section 
show how this precise topological field theory arises as the ground state of a specific Hamiltonian. We 
note in addition that if one specializes these results to an a belian ga uge field, the topological field 
theory obtain ed i s known as a (two-dimensional) BF theory 109lll24l ]. Recently, Freedman, Nayak 



and Shtengel |l25l | have discussed (reasonably local) lattice models of interacting fermions and bosons, 
which they argue have topological ground states with some of the topological algebraic structure of the 
doubled non-abelian Chern-Simons theories. It is reasonable to expect that the universal long-distance 
structure of the wave functions of these topological ground states, at least deep in the topologically- 
ordered phase, has the same structure of the wave functions we discussed in this section. 

Since we have an explicit Hamiltonian Eq. ()5.26() (plus the analogous term for the B gauge field), 
we can go beyond the topological field theory describing the ground state. Because of the commutation 
relations ()5.25|) . E'^ acts like a creation operator. However, in a non-abelian gauge theory, E and E'^ 
are not gauge invariant, s o th e appropriate states are slightly more complicated than just 
The simplest candidate is [llOl ] 

i)''{x,A,B) = E^\x) j[Dg]Tr T''g{x)T''g-\x) e'^^^A^B) _ (5 35) 

Once this amplitude is squared and B is integrated over, one obtains a gauge-invariant probability. In 
the 2 + 1-dimensional picture, one can think of one of these states as a Polyakov loop, which intersects 
two-dimensional space at a single point x. The commutation relation Eq. (|5.25l) shows that this is 
indeed an eigenstate with a gap proportional to k. A more thorough treatment, taking into account 
the measure of the path integral, shows that the gap is shifted to e^(fc + 2c4)/(47r), where ca is the 
quadratic Casimir of the adjoint representation of the Lie algebra of G For our purposes, the 

important point is that there is indeed a gap. 

We have thus seen that correlators of Wilson loops in the ground state of strongly-coupled Yang- 
Mills theory are given by a topological field theory. We have also seen that the theory has a gap, so it 
is indeed in a topological phase. The last thing we would like to discuss is the ground-state correlators 
for operators not in the topological theory. In other words, we wish to consider the full set of physical 
operators (i.e. those annihilated by Qbrst), including those which are given by Qbrst acting on 
something else. Because our Hamiltonian involves only the electric and not the magnetic field, it 
is not Lorentz invariant. Thus the existence of a gap does not immediately require that correlators 
decay exponentially. However, in this theory, they do, as implied by the results of jll9l |. We noted 
above that proper quantization of gauged WZW models requires introducing fermionic ghosts, and a 
fermionic operator Qbrst- A fermionic symmetry suggests the appearance of supersymmetry, and 
indeed the G/G topological field theory can also be obtained from a supersymmetric field theory where 
the supersymmetry charge is "twisted" into Qbrst jl09' 126T 271 128| |. The supersymmetric field theory 



also describes the correlators of the full theory, not just the topological subsector. For the case where 
G = SU{N) at level k, the appropriate two-dimensional field theory is a (twisted) supersymmetric 
sigma mod el, where the bosonic fields take values on the "Grassmannian" manifold U{N + k)/[U{N) x 
U{k)) |ll9l ]. It is likely that the other simple Lie groups end up giving supersymmetric sigma models 
on the analogous Grassmannians. 
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This means we have now come full circle! Wc started without gauge fields, and found that the 
equal-time correlators are those of the principal chiral model, a two-dimensional sigma model with a 
curved target space. These correlators are exponentially decaying. We then introduced gauge fields, in 
the hope of finding a quantum critical point. After this lengthy discussion, we have ended up showing 
the ground state is described by a supcrsymmctric two-dimensional sigma model with a curved target 
space, the Grassmannian. This means that the correlators here are exponentially decaying as well, 
and the theory is in a topological phase whose properties are encoded in the wave function x[A,B]. 
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APPENDIX A: OPERATORS OF THE QUANTUM LIFSHITZ FIELD THEORY 

In addition to the products of field operators </?(x), in what follows we will be interested in two 
types of local operators: charge and vortex operators. The charge operators are 

Onix) = e-^" ^(^) (Al) 

where n G Z. This operator creates a boson coherent-state which we will refer to as a charge n 
excitation. The vortex operators are 

i I (fz ifiz) U(z) 
Om{S) = eJ (A2) 

where 

^{z) = marg(z — x) (A-3) 

where < arg(z — x) < 27r is the argument of the vector z — x (with a branch cut defined arbitrarily 
along the negative x axis). The action of the operator Om{x) on an eigenstate of the field operator 
I [(/?]) is simply a shift 

i I (fz ip(z) U(z) 

eJ M) = Mx)-vm (A4) 

In other words, it amounts to s singular gauge transformation. Therefore, its action is equivalent to 
coupling the field (/? to a vector potential whose space components A satisfy 

dz ■ A[z\ = 27rm (A5) 

7 



38 



for all closed paths 7 which have the point x in their interior, and zero otherwise. In particular, the 
wave function of the state resulting from the action of the vortex operator on the ground state is: 



^m[x] = {[ip]\d^mix)\vac) 



1 



A 2 



(A6) 



where A is any vector field which satisfies Eq. (jASf) . The (equal-time) ground state expectation value 
of a product of vortex operators with magnetic charges {mi}, i.e. the overlap of the state with k 
vortices at locations x/ and magnetic charge m/ with the vortex-free ground state wave function, is 
therefore 



(vac|Omi(^i) • • • Omfc(^fc)|vac) = 4 I 'Dip e 



K d^z [Vip - A 



where Z is given by Eq. (|2.12j) . The vector potential in Eq. ()A7|) satisfies 



z-xi) 



(A7) 



(AS) 



1=1 



This result is equivalent to the expectation value of the vortex operator in the 2D classical c = 1 
compactified free bose field discussed extensively by Kadanoff jo^ (see also [48.]). 
The boson propagator of this theory, in imaginary time t, is 



G{x - x,t -t') = {ip{x,t) ip{x,t')) = / — 



duj 



d'^q e 



iuj{t — t') — iq- {x — x) 



2tt J (27r)2 



+ K? {q"^Y 



(A9) 



which has a short-distance logarithmic divergence. From now on we will use instead the regularized 
(subtracted) propagator 



Grcg(x,t) = G{x,t) - G{a,0) 



8ttk 



X P \ / IxP 



where a is a short-distance cutoff and r(0, z) is the incomplete Gamma function 



r(o,z) 



The regularized propagator has the asymptotic behaviors 



(AlO) 



(All) 



Greg(x,t) = < 



where ln7 = C = 0.577 ... is the Euler constant. 

The correlation functions of the charge operators are 




(A12) 



{On{x, t)^On{x', t')) = e '^-g(^ -X',t-t') 



(A13) 



At equal (imaginary) times, |t — t'| ^0, it behaves like 

(O„(x,0)tO„(f' ,0)) = 



X — X 



n 

Attk 



(A14) 
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which imphes that the operator On has (spacial) scahng dimension 

^2 



n 



Sttk 

For |x — x'l — > a, its asymptotic behavior is instead given by 



(A15) 



(O„(0, t)tO„(0, t')) = (4^^^) (A16) 

This behavior is manifestly consistent with a dynamical critical exponent z = 2. 

It is straightforward to show by an explicit calculation of the overlap of Eq. (jASI) , which is completely 
analogous to the classical vortex correlation functions of the 2D Gaussian model (69|], that the (spacial) 
scaling dimension of the vortex of magnetic charge m is, 

Am = 27rKm^ (A17) 



APPENDIX B: C/(l) GAUGE-THEORY FOR THE QUANTUM SIX- VERTEX MODEL 

In this Appendix, we describe the quantum six-vertex model in the language of gauge theory. The 
quantum eight- vertex model will be described in appendix [O We will follow closely the gauge theory 
description of the quantum dimer model, which is described in detail in ^^l- We note that this gauge 
theory is not simply an abelian version of the gauge theory discussed in section IVl 

We define link variables Ei{x), where x labels the vertices and i = 1, 2 indicates the direction; the 
unit vector in direction i is denoted by ej. The link variables are integer valued, and can be viewed 
as the eigenvalues of angular momentum operators which we will also denote by Ei{x). We assign the 
values El = 1 and £^2 = 1 to the right and up going arrows respectively, while for the left and down 
going arrows we have Ei = —1 and E2 = —1. An example is given in figure [Tj In the Hamiltonian, 





E2{x) = 1 


Ei[x-ei) = 1 


X 


Ei{x) = -l 


E2{x-e2) = ~l 





FIG. 7: The 'electric' fields associated to a particular vertex 



we need a term which will lead to a restriction to those states in which all the £'i(x)'s have values ±1. 
Such a term is 

in the limit /ci ^ 0. 

The six-vertex model also requires that the same number of arrows point in and out at a vertex. 
This is precisely a lattice version of Gauss' law. Defining the lattice divergence as A'-E-lx) = Ei{x) — 
Ei{x — ei) + E2{x) — £'2 (x — 62), one easily finds that only for the six vertices of type a, h and c, we 
have 



(B2) 
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The lattice differentiation is defined by = f{x) — f{x — e,). In tlie following, when dealing 

with sums over the plaquettes of the lattice, we will frequently use Aj/(x) = A't/(x) = /(x + Cj) — /(x). 
Of course, the constraint has to commute with all the terms in the Hamiltonian. In the following, we 
will find that this is indeed the case. 

The two main ingredients in the gauge theory description arc the flip term which flips the flippable 
plaquettes (see below) and a potential term, which give a finite weight to (only) the fiippable plaquettes. 
Flippable plaquettes are those which have both = = 0. Pictorially, a flippable plaquette here is 
one where the arrows around the plaquette point either all clockwise or all counterclockwise: 



-r^ -t-r . (B3) 

In terms of the electric field, this can be written as 

= ^ Y,i^2Eif{^iE^f{E^ - E^f , (B4) 
□ 

where the sum is over all the plaquettes of the lattice. The factors (A2-Bi)^ and make sure 

that the arrows on opposite links of the plaquette arc anti-parallel. The factor (^Ei — £^2)^ checks if 
two arrows on one vertex are both pointing clockwise or both counterclockwise. 

As the £^i(x)'s have the integers as their eigenvalues, the canonically conjugate operators, ai{x) are 
phases, i.e. < ai{x) < 27r. Using the commutation relations 

[aj{x),Ej,{x )] =iSjf5g^g> , (B5) 

it is easy to show that the operators e^^^ C^) act as raising and lowering operators on the £"s, and thus 
e.g. the operator e"^*"^^^) will fiip the arrow pointing outward from x to the arrow pointing inward. 
We can use these raising and lowering operators to write the term in the Hamiltonian which flips the 
flippable plaquettes 

Ht = -2t J2 cos(2 Eaajix)) , (B6) 
□ 

where Yla'^ji^) = •^102 (^) — A2ai(x) is the oriented sum of the a's around a plaquette. The total 
Hamiltonian of the gauge theory version of the 6-vertex model is therefore 

He^ = HE + Ht + Hy . (B7) 

As usual, the Rokhsar-Kivelson point is located at t = V. 

We will proceed by going to the dual formulation of the theory, and show that the theory is 
equivalent to a height model (which is well know). Doing the duality basically amounts to solving 
the (electrostatic) constraint. In the process, it gets replaced by a magnetic constraint. To solve the 
constraint, we introduce the new variables S{r), which live on the sites of the dual lattice (or plaquettes 
of the direct lattice); these operators have the integers as their spectrum. In addition, we need the 
fields Bi{f), which live on the links of the dual lattice. We can now write the "electric" fields Ei as 
follows: 

Eiix) = eij[A-Sir) + Bj{r)) . (B8) 

Substituting this in the constraint ^^E- = gives the 'magnetic' constraint eij/^-B^(f) = 0, so is 
curl free, and can be written as a gradient. But, as the there are no sources, we can do even better, as 
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becomes clear when we interpret S{r) as a height variable which lives on the plaquettes of the direct 
lattice. 

We will first recall the known fact that the configurations of the 6-vertex model can be mapped 
onto height configurations. The rules are as follows. First, pick a reference site, and give it a reference 
height, say S{0) = 0. Then, if one crosses an outgoing arrow clockwise (both seen from the vertex), the 
height decreases by one, while crossing an incoming arrow (again in a clockwise manner), the height 
increases by one. As all the vertices have two incoming and two outgoing arrows, this indeed gives a 
consistent height configuration. So, as an example, around an a vertex we have 



1 








-1 



We now assume that the S'(r)'s appearing in Eq. (|B8|) can in fact be interpreted as the heights of the 
plaquettes. Because of the duality Eq. (|B8|) . the Bi are now completely determined by the Ei, because 
they determine S via the height rule. So we found that interpreting the S'(r) as heights amounts to 
picking a gauge for the Bi (r ) . Combining the height rule of the previous paragraph with Eq. (jBSf) , we 
easily find that Bj{r) = 0. 

To complete the duality transformation, we need to transform the flip term. This will involve 
the canonically conjugate variable to S{f). Let us call this the momentum P{f), which satisfies 
[P{r), S{r')] = i6ff^f'. Again, acting with e^^^^ on the plaquette at r will increase the eigenvalue of 
S{r) by one. Flipping a plaquette changes S by ±2, so we find that we can write the flip term of 
the Hamiltonian as Ht = — 2t ^^cos(2P(r)). In other words, we find that the circulation around the 
plaquette at r is given by Yla'^ji^) — P{^- 

In principle, we could go on to describe the eight-vertex model in a similar fashion. The main 
difference with the six-vertex model is that the constraint HB2() is no longer satisfied. One can introduce 
a matter field which has no-zero values on the d vertices. In addition, a flip term which flips every 
plaquette has to be constructed. This flip term has to commute with the new constraint. It turns out 
that this is indeed possible, but the flip term will involve the conjugate of the matter field. Also, one 
would need a potential term which gives weights plaquettes according to the vertices present on the 
plaquette. However, as there is a more natural gauge description, which is based on the Z2 symmetry 
of the eight-vertex model, we will discuss and use that description of the quantum eight vertex model 
in the next appendix. 

APPENDIX C: A Z2 GAUGE THEORY FOR THE QUANTUM EIGHT VERTEX MODEL 

In this appendix, we will discuss a Z2 gauge theory which can be viewed as an extension of the 
Kitaev model which incorporates vertex weights differing from unity. The model we will discuss is of the 
Rokhsar-Kivelson type, but it is not the simplest quantum generalization of the classical eight-vertex 
model, as we pointed out in section IIVI 

1. The Z2 gauge theory 

In this model, spins living on the bonds of the square lattice are the degrees of freedom. Thus, on 
every link of the square lattice, we define a Pauli algebra of 2 x 2 Hermitian matrices cr"(x), where 
a = 1, 2, 3 labels the three Pauli matrices, and for a lattice site x, we denoted the orientation of the link 
by j = 1,2 (1 =horizontal and 2 =vertical). (Thus, the degrees of freedom live half-way between the 
lattice sites x and x + ej, where ej is a unit vector along the direction j.) In what follows we will take 
the states | t ) (an "up" spin) and | i ) ( a "down" spin) as the states which diagonalize (instead 
of o"^, as it is customary). The relation with the eight-vertex model is simple. Up-spins correspond 
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Gl[x-er] 




o<i{x) 



FIG. 8: Notation for the c's; a = 1,2,3 is the Pauh matrix label. 



to arrows pointing up or to the right, while down spins correspond to arrows pointing down or to the 
left. Around each vertex, the number of up-spins has to be even. In this section, we will denote the 
vertex by x and the associated plaquette by its south-west corner x. 

For details, see Ref. ^31-5 7.] . The constraint can written in a simple form^'' 

a\{x) a\(x — e\) a\{x) a\(x — e'2) = \^ Vx (CI) 

The term which flips all the arrows around a plaquette can be written in terms of u^'s 

ii^fiip = -^o\{x) fT?(f +62) al{x) cTi(x + ei) (C2) 

This flip operator commutes with the constraint, because a vertex and a plaquette have either or 2 
common bonds. Hence, both operators can be diagonalized simultaneously. In the model considered 
by Kitaev [i^j, equal weight is given to all types of vertices, so no term in the Hamiltonian is required. 
To go beyond the point a = 6 = c = (i=l, we need a term which will weight the plaquettes according 
to the types of vertices present. Thus, we will need to introduce operators which can discern among 
the various vertices. In addition, the "weight" of the vertex depends on its position in the plaquette 
under consideration. The notation that we use is summarized in figure |S1 

Let us start by giving the terms which give a non-zero contribution for a, 6, c and d vertices respec- 
tively (but zero otherwise). Let us define the "vertex magnetizations" 

6*1 ) +(T2(X -62)) 

ei) - a\{x-e2)) 
ei) + <T2(x - 62)) 

ex)-a\{x-e2)) (C3) 

With this notation, the projectors onto the vertices a, 6, c and d are just the squares of the vertex 
magnetizations of Eq. ()C3|) : 

Va = S^, V, = S,\ Vc = S,^ Vd = Si (C4) 



Sa{x) = ^{a\{x) + al{x) + a\{x - 

Sb{x) = ^{a\{x) - alix) + a\{x - 

Sc{x) = ^ {(t\{x) - al{x) - a\{x- 

Sd{x) = ^{(j\{x) + a\{x) - a\{x - 



Throughout we use the upper label to indicate the Pauli matrix and the lower label to indicate direction. 
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It is straightforward to show that these operators act as projection operators on vertices of type o, b, 
c and d respectively, i.e. they yield 1 when acting on the corresponding vertex and zero otherwise. It 
is elementary to check that 

ra + rb + Vc + Vd = i (C5) 

where I is the identity operator. 

We can now write down the potential term which assigns weights to the plaquettes in the same way 
as done by the projectors of section IIVI Flipping a plaquette will change a d vertex to an a vertex if 
the d vertex is in the south-west or north-east "even" corner of the plaquette. A d at the other "odd" 
corners will go to a 6 under the flip. In short, dg <-> ag, do <-> bo- For c's this is opposite, namely 
Ce be, Co do- The potential term will involve all four vertices around a plaquette, and thus we 
need to distinguish between the different position in the plaquette. This is achieved by the following, 
albeit rather cumbersome, term 



= V (-Vaix) + -Mx) + -Vc{x) + ydx)) 

^ \a b c d I 



-Va{x + ei) + ^Vb{x + ei) + -Vc{x + d) + ^Vd{x + ei)] x 

a bed) (C6) 

-Va{x + e2) + ■T'Pb{x + e2) + -Vc{x + e2) + iPd{x + e2)] x 
a b c d J 

-Vaix + ei + 62) + jVh{x + ei + 62) + -Vc{x + ei + 62) + '^Vd{x + ei + €2] 
abed 



The potential term HC6|) assigns potential energies to the plaquettes in the same way as is done by the 
projectors of section HVl That is, the potential is the product of vertex weights obtained by flipping 
the plaquette, divided by the product of the vertex weights of the plaquette itself. As mentioned 
before, the two-body terms only couple spins on the same sublattice. The most non-local terms in 
the potential energy term consist of eight-body interactions. Of course, a plaquette potential energy 
term of this sort is needed, if one assigns weights to plaquettes^^. Note that for a = b = 1 and c = d 
(or e = 2)1 the eight-body terms cancel each other. However, there will remain four and six-body 
interactions, which will mix different sublattices. 

The total Hamiltonian for the Z2 gauge theory of the quantum eight vertex model is 

HqSv = Hy + Hf^ip (C7) 



where the two terms are given by HC2p and HC6|) . The states in this model have to satisfy the constraint 
()C1|) . We can check these formulas in a few limits. For a = b = e = d= 1, each factor in (|C6|) will 
be 1 when acting on states satisfying the constraint. Hence, the potential term became the identity 
operator, which merely results in an energy shift, and thus we find back the Kitaev model, as we must. 

Another interesting limit is d — > 0. As was discussed in section IIV Al this limit gives a slight 
generalization of the six-vertex model, in the sense that plaquettes which have = = 1 or 
nd = nd = 2 will have finite energy, and they can be interpreted as static defects. As can be seen 
from the potential term lC6l plaquettes which have > will be suppressed as they receive infinite 
energy. The flip term has to be modified, because we need to have a flip term which commutes with 
the constraint, which has become a stronger statement in the six- vertex case, namely (c.f. Eq. HB2|) ) 

al{x) - a\{x - ei) + al{x) - al{x - 62) = Vf . (C8) 



In spite of the appearances, Hv as given in Eq. IIC6t respects rotational invariance; the apparent asymmetry is due to 
the use of the vertex weights as labels. 
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The flip term which preserves the six-vertex constraints is 

^flip,6v = -y^('7r(^) cr^jx + ei) af{x + e2) o-^(x) + h.c.) , (C9) 

X 

where the raising and lowering operators (in the representation we use) are given by 

a^ = ^{a^^ia^) (CIO) 
To make contact with the flip term for the eight-vertex model, we rewrite Eq. H4.12() as 

-H"flip,q8v = - yZi^i(^) + ^r(^)) (^^2^ + ei) + {x + ei)) X 

i _ _ (Cll) 

[aXix + 62) (Jx (a^+ 62)) ('7^(2;) + (^)) • 

The flip term for the six-vertex model is therefore precisely the flip term for the eight-vertex model 
minus the terms which cause the six-vertex constraint to be violated. It is easily checked that the 
six- vertex flip term Eq. HC9|) commutes with the constraint ()C8|) . We thus find that on the level of 
the wave function, the limit d ^ is smooth, as the amplitude of the configurations which contain d 
vertices goes to zero. In addition, for d 7^ 0, the flip term commutes with constraint HC1|) . while for 
(i = 0, the flip term commutes with the ^7(1) constraint HC8|) . Thus, the symmetry is enhanced from 
Z2 for d 7^ to f/(l) for d = 0, as was to be expected. 



2. The dual of the gauge theory 



r — ei 
• 


r 

• 

X 


\ < 

f - 6*1 - 6*2 
• 


\ \ 

f - 6*2 
• 



FIG. 9: The dual lattice sites are labeled by r. 



We now have a Rokhsar-Kivelson generalization of the eight-vertex model, in a gauge-theory lan- 
guage. We can use this representation of our model to study the various phases. However, this is more 
easily done in a dualized version, as the dual takes the form of an Ising model. In the dual picture, the 
spin degrees of freedom will live on the sites of dual square lattice, i.e. the centers of the plaquettes 
of the direct lattice. Thus, we will label by r the site of the dual lattice on the center of the plaquette 
labeled by x (its SW corner). Of course, the potential term in the dual language will still be quite 
formidable. We will denote the dual Pauli operators by and r^. To start with the flip term, the 
product of cr^'s around a plaquette becomes on the plaquette [sillsTt 

T^(f) = al{x) al{x + ei) a?(f +6*2) al{x) (C12) 

To see what happens with the constraint and the projector operators defined by Eq. HC3|1 and Eq. 
()C4() we need the dual form of the 's living on the links. In term of the dual variables r^, and using 
the notation of Fig. |3 the 's are given by 

a\{x) = T^{r)T^{r-e2), ct^ (x) = (f ) (f - ei ) (C13) 
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We thus easily find that the constraint is automatically satisfied (again, going to the dual picture 
amounts to solving the constraint). Also, it is trivial to show [sills^ that the inverse relation, i.e. to 
express the dual lattice operators in terms of the operators of the original lattice is 

r'ir) = n ^\^) (C14) 

where {£} is a set of links of the direct lattice pierced by a path r(r) on the dual lattice ending at the 
dual site r (but which is otherwise arbitrary); see Fig. 1101 
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FIG. 10: The dual path T. 



Note that we need to choose the spin on one of the plaquettes; all the others are subsequently 
determined by the cr^'s on the links. In terms of the dual variables, the projection operators for site 
X, defined by Eq. HC4|) . take the form 



Vaix) = ^{1+A{f^+ B{r) + C(f)) 
V,{x) = \{l-A{f^-B{r)+C{r)) 
V^{x) = ^ (1 - Air) + B{r) - C(f)) 
V,{x) = ^{l+A{f^-B{r)-C{r)) 



(C15) 



where A, B and C are given by 



A{r) = T^{f- ei) T^(f - 62) 

B{r) = T^{r)T'^if- 61-62) (C16) 
C(r) = T^(f) r^(f - ei) r^(f - ei) r^(f - ei - 62) . 

The interaction in terms of these projectors has the same structure as in the spin representation of 
the classical eight-vertex model: it consists of two-body terms on interpenetrating sublattices, and a 
four-body term coupling the two sublattices. All the two-body interaction terms will only couple spins 
on the same sublattice. The four- and six-body interaction terms (of which there are many!), couple 
the sublattices. The same holds for the eight-body term, naturally. The dual form of the theory is 

^q8v,dual = ^V,dual " X] '^^ ' ('^^ 

r 

where ^^v,duai is given by (|C6() . but now with the projectors given in Eq. (|C15() . Thus, formally this 
theory takes the form of a (multi-spin) Ising model in a transverse field. However, the two-body 
interactions only couple spins on the same sublattices, together with the multi-spin terms conspire to 



46 



change the quantum critical behavior from the conventional z = 1 Lorentz-invariant criticality of the 
standard Ising model in a transverse field to the z = 2 quantum critical behavior discussed in the rest 
of this paper. 

Now that we found the dual version of our gauge theory, we would like to discuss the limits 
a = b = c = d = l and d = 0. Again, the first limit brings us back to the Kitaev point, because 
the potential term becomes the identity operator again, and we are left with the very simple spin flip 
term of Eq. HC17|) . Hf = — Ylr'^^i^)- The limit d — > is however more complicated in this dual gauge 
theory. First of all, we now do need a constraint, which was not present for d ^ 0. Moreover, the flip 
term now only can act, depending on the surrounding spins. Let us start by dualizing the constraint 
Eq. ()C8() . which results in 

(r3(f)-r3(f-ei-e2))(r3(f-ei) + r3(f-e2)) =0 Vf . (CIS) 

Obviously, the eight-vertex flip term r^(r) does not commute with this constraint. To find a flip term 
which does commute with the constraint, we dualize the six- vertex flip term <\(J9\i . which results in 

^^flip,q6v = ^Hr) (1 - r\r + ei)T^{f+ ea)) (l + T^{f- ei)r3(f + es)) (l + T^{f+ ei)r=^(f - es)) . 

(C19) 

The factors in bracket can be seen to give a non zero result only on plaquettes which are flippable. 
Hence, this flip term commutes with the constraint (|C18() . Of course, this can also be checked explicitly. 
Apart from a factor r^, there are factors depending on the r^'s coming form both the flip term and 
the constraint. The signs in this product conspire in such a way to render the commutator zero. 
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